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The  ionosphere  is  the  region  of  the  atmosphere  in 
which  sufficient  numbers  of  free  electrons  and  ions 
exist  to  affect  radio  wave  propagation.  This 
ionization  is  mostly  caused  by  solar  radiation 
interacting  with  the  Earth's  atmosphere.  The  electron 
density  distribution  of  the  local  ionosphere  changes 
with  solar  variability,  such  as  day  and  nighttime 
difference,  seasonal  variations,  sunspot  number,  etc. 

Although  free  electrons  and  ions  are  distributed 
randomly  on  small  spatial  scales,  they  often  act  upon 
siectromagneti c  waves  as  a  continuous  refracting  medium 
because  of  their  statistical  uniformity.  In  the  high 
frequency  (HF)  range  the  ionosphere  may  be  considered 
as  composed  cf  cniy  free  electrons,  since  the  ions  have 
heavy  mass  and  are  relatively  immobile. 

The  collisions  free  elections  make  with  other 
particles  such  as  ions,  molecules  and  atoms  have  a 
viscous  damping  effect  on  radio  waves  which  is 
equivalent  to  ohmic  heating  in  a  conductive  medium. 

This  is  the  major  reason  fcr  radio  wave  absorption  in 


one  f "nos p here. 


» 


The  Earth's  magnetic  field  complicates  the 

refracting  properties  of  the  ionosphere.  The  presence 

of  a  constant  magnetic  field  makes  the  medium 

anisotropic  and  complicates  tne  relationship  between 

the  electric  polarization  P  and  electric  field 
— 

'intensity  E, 

P-  e a  M  E  (1.1) 

where  M  is  the  3x3  susceptibility  tensor  of  the 
ionized  medium  and  e0  is  the  free  space  permittivity. 

As  an  anisotropic  medium,  the  ionosphere  exhibits 
many  interesting  properties.  For  special  cases, 
complete  penetration  of  the  layers  which  would 
generally  be  reflecting  is  possible  through  ordinary  tc 
extraordinary  mode  wave  coupling.  This  occurs  when 
pro  pagatim  is  near  the  critical  ar.gia  of  incidence 
which  lies  in  the  magnetic  meridian.  Total  energy 
deposition  of  the  coupled  wave  may  occur  near  the 
resonance  region  using  this  mechanism.  In  the  f  it  are 
this  may  prove  to  be  an  efficient  heating  mechanism  fcr 
ionospheric  modification  experiments  as  well  as  a 
viable  explanation  fcr  many  experimental  observations  . 
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A  theoretical  and  computational  examination  of  this 
phenomenon  is  made  In  this  study. 

1.2  -  IONOSPHERIC  HEATING  EXPERIMENTS 

High  power  HF  heating  experiments  began  in  the 
1970s  and  produced  many  interesting  and  often 
unexpected  results  [1,2].  The  motivation  of  the  early 
experiments  was  to  understand  the  nature  of  the 
ionospheric  plasma.  Experimental  observations  have 
been  summarized  in  numerous  references  [ 3 , 5 , 5 , 6 , 7 , 8 , 9  , 
10,11,12,13].  Among  the  most  interesting  observations 
derived  from  these  experiments  are: 

-  Plasma  Instabilities  associated  with  the  non¬ 
linearity  behavior  of  the  medium. 

-  Generation  of  short  scale  i  rr  egul  ar  i  t  i  es  . 

-  Generation  of  large  scale  irregularities. 

-  Airgiow  excitation. 

-  Nonlinear  demodulation  (ELF,  VLF  generation). 

-  Production  of  extrathermal  electrons. 

These  experimental  findings  may  have  their 
greatest  impact  on  long-range  communications  systems. 
For  example,  scattering  from  sm3li-scaj.e  irregularities 


(0.5  -  5m  size)  have  been  observed  which  may  be  useful 
for  establishing  long  range  VHF  and  UHF  communication 
links  [33,34].  Large  scale  Irregularities  were  found 
to  scatter  HF  signals  [14,15,16,17]  which  might  be 
useful  for  injecting  HF  waves  into  ionospheric 
waveguides  for  so-called  duct"  mode  propagation  [18,19, 
20,21,22].  Nonlinear  demodulation  capability  of  the 
ionosphere  has  also  been  reported  [23,24].  When  the 
heater  wave  is  amplitude  modulated  with  a  frequency 
between  0.5~10KHz,  the  heated  volume  irradiates. waves 
at  the  modulating  frequency  which  may  be  received  at 
far  distances.  This  is  believed  to  occur  because  of 
the  modulation  of  the  natural  ionospheric  currents 
(dynamo,  polar  electrojet,  equatorial  electrojet,  etc.) 
passing  through  the  periodically  varying  heated  region. 

The  potential  of  HF  ionospheric  heating/ 
modification  has  been  well  recognized.  However,  tc 
date,  most  of  the  experiments  have  been  performed  with 
vertically  directed  antennas.  Both  theoretically  and 
experimentally,  very  little  attention  has  been  paid  to 
oblique  propagation.  This  work  shows  that,  for  certain 
angles  of  incidence,  oblique  propagation  should  provide 
an  efficient  way  of  depositing  energy  into  the 
ionospheric  piasma.  This,  then,  should  enhance  the 


ut.llity  of  HF  radiowave  heating  for  communications  and 


basic  plasma  research  purposes. 

1.3-  RESEARCH  APPROACH 

Electromagnetic  wave  propagation  in  a  plane 
stratified  ionospheric  plasma  and  its  possible 
application  to  high  power,  high  frequency  (HF)  heating 
experiments  are  investigated.  Theory  shows  that  mode 
coupling  from  ordinary  (o- )  to  e xtraor di nar y  (x-)  waves 
in  the  vicinity  of  the  plasma  resonance  level  is 
possible  if  certain  conditions  are  met.  This  may 
result  in  deposition  of  radiowave  energy  into  the 
ionospheric  plasma.  The  conditions  under  which 
coupling  is  most  effective,  and  the  amount  of  coupling, 
are  investigated  for  several  ionospheric  models  using 
the  numerical  full  wave  solution  technique  and  the 
roots  of  the  3ooker  quartic. 

The  theory  of  wave  propagation  in  the  ionospheric 
plasma,  the  Booker  quartic  and  the  full  wave  solution 
technique  are  presented  in  Chapter  2.  Computational 
results  which  demonstrate  the  theory  and  coupling 
efficiencies  for  several  directions  of  propagation  are 
presented  In  Chapter  3.  Summary  of  results, 
conclusions  and  recommendations  are  presented  in 


1 


1 
a 


$ 

*  *1 


-5- 


Chapter  4.  This  study  cieariy  shows  that  coupling  from 
o-  to  x-waves  may  result  in  significant  absorption  of 
radiowave  energy  in  a  very  narrow  region  just  below  the 


level  at  which  the  wave  frequency  equals  the  plasma 
frequency.  This  heating  mechanism  should  have 
important  implications  for  using  ionospheric  high-power 
radiowaves  to  produce  l rregul ar i t i es  from  which 
scattering  of  other  radiowaves  for  communications 
purposes  may  be  possible. 


CHAPTER  2 


REVIEW  OF  RELATED  THEORY 

In  this  chapter  the  theory  of  HF  electromagnetic 
plane  wave  propagation  in  a  plane  stratified 
ionospheric  plasma  is  reviewed  [25,26,27,28,29,30.]. 

The  theory  provides  the  basis  for  understanding  HF  wave 
propagation  phenomena  through  related  numerical  full 
wave  sol ut i ons . 

2.1  -  ASSUMPTIONS 

-  The  electromagnetic  fields  are  assumed  to  arise 
from  harmonic  plane  waves.  The  time  dependence 
exp(lwt)  is  suppressed  for  all  field  variables. 
The  rationalized  M.K.S.  units  are  used. 

The  ionosphere  is  assumed  to  be  horizontally 
stratified. 

As  a  magnetoionic  medium,  the  ionosphere  is 
assumed  to  consist  of  free  electrons  which  are 
distributed  locally  with  statistical  uniformity. 
The  effects  of  heavy  ions  are  neglected. 

Cold  plasma  conditions  are  assumed  to  exist, 
and  the  possibility  of  energy  loss  caused  by 
generation  of  other  types  of  waves  such  as 
plasma  waves  and  acoustic  waves  is  ignored. 


Definition  of  the  symbols  and  basic  quantities 


i  are  listed  in  Table  2.1. 

! 

2.2  -  BASIC  EQUATIONS 

The  electromagnetic  fields  in  the  ionosphere  are 
governed  by  the  Maxwell's  equations: 


V  •  D  -  0 


(2.1) 


V  x  E  -  1  i  u>  p  0  H 


(2.3) 


V  •  B  -  0 


(2.2)  7  x  H  -  i  u  D 


(2.4) 


If  the  z-axis  is  taken  to  be  the  wave  normal,  Eqs .  2.3 
2.4  Become 


3_ 

3z 


a)  Ho  H  x 


(D  D  y 


3_ 

3  z 


Ex 


-  i  ID  D0  H  y 


-  i  a  Dx 


(2.5 


H  z  =  0 


Dz  »  0  leads  to  Ez  »  0  in  an  Isotropic  medium.  How¬ 
ever,  in  an  anisotropic  medium  such  as  the  ionosphere 
this  is  not  true,  and  in  general  the  longitudinal 


component  Ez  is  non-zero. 


Tab le  2.1  -  Symbols  and  Notation 

velocity  of  light  in  free  space 
free  space  permittivity 
free  space  permeability 

/^T 

electron  charge 
electron  mass 
wave  number 
Boltzmann  constant 
electron  temperature 
refractive  index 
electron  density 
electric  field 
electric  polarization 
electric  displacement 
magnet i c  field 
magnetic  induction 
n0  H 

wave  frequency 
plasma  frequency 


gyro  frequency 


Earth's  magnetic  induction 


plasma  frequency  parameter 
normalized  gyro  frequency  parameter 
normalized  collision  frequency 
3x3  susceptibility  tensor 
x,y,z  direction  cosines  of  Y 
1  -  iZ 

collision  frequency 

collision  frequency  of  monoener get i c  electrons 

effective  collision  frequency 

complex  collision  frequency 

Dingle  integral  of  order  n 

wave  polarization 

transverse  component  of  Y 

longitudinal  component  of  Y 

=*  1-Y2  resonance  condition 

angle  of  incidence 

angle  of  incidence  in  free  space 

magnetic  dip 

intrinsic  impedance  of  free  space 
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H  is  often  used  as  a  measure  of  magnetic  field 
intensity  and  is  defined  as: 


H  =  n«  H 


(2.6) 


where  n0  is  the  intrinsic  impedance  of  free  space. 

H  has  the  units  of  electric  field  intensity  which  makes 
it  more  convenient  for  this  application.  Using  this 
definition  Eqs .  2.3,  2.4  become 


V  x  E  =  -  i  k  h 


-*  i  k 
V  x  H  =  D 
£  « 


(2.7) 


(2.8) 


In  the  ionospheric  plasma  the  equation  of  electron 


motion,  arising  from  excitation  of  an  HF  electro¬ 
magnetic  wave,  defines  the  polarisation  state  of  the 
plasma  according  to  Eq.  2.9  [26] 


-e0  XE=P(l-iZ)+iPxY 


(2.9) 


where  X,  Y,  Z  are  defined  in  terms  of  basic  quantities 


IJ.* 


VXYJW 


(2.10) 


e  B  o  =  oiH 


(2.11) 


2  =  - 
0) 


(2.12) 


The  variable  X  is  proportional  to  the  electron  density. 
Plasma  frequency  oi^  Is  a  characteristic  frequency  of 
the  plasma  at  which  electrons  tend  to  oscillate.  The 
variable  Y  is  the  normalized  gyrof requency ,  and  its 
direction  is  ant i -> par  ail  el  to  the  superimposed  magnetic 
field  $0 .  The  gyr of r e q uen c y  is  another 
characteristic  frequency  of  the  magnetoionic  plasma  and 


is  the  frequency  at  which  free  electrons  oscillate  in  a 
circular  motion  about  the  magnetic  field  Bg .  The 
variable  Z  is  the  normalized  collision  frequency. 

If  (l,  m,  n)  are  the  direction  cosines  of  Y, 

E q .  2.9  becomes 


i  n  Y 


i  m  Y 


“  e  o  X  Ey 


-i  n  Y 


i  l  Y  Py  (2.13) 


i  m  Y  -i  l  Y 


Inverting  this  3x3  matrix  gives  the  classical  form  cf 
the  constitutive  relations  r2Bl 


Uses  .of  the  parameter  Z  suggests  that  there  is  a 
single  effective  collision  frequency  which  can  be 
substituted  in  the  constitutive  relations.  Actually, 
because’  of  the  dependence  of  the  collision  frequency  on 
thermal  energy  of  the  electrons  there  is  no  single 
collision  frequency  applicable,  except  under  certain 
conditions.  Sen  and  Wyller  (i960)  [31]  developed  a 
modified  form  of  the  constitutive  relations  for  weakly 
ionized  gases  which  accounts  for  the  thermal  energies 
of  electrons  using  Boltzmann’s  transport  equation.  As 
a  result  one  obtains  three  different  collision 
frequencies  gn  (n-1,2,3)  [32]. 

8n  =  un  if  CX  ( ^n)  1  -  i  u)n  (2.15) 

*  7  vm  ~  vm 

where 

w  i  *  w  <jj  2  =*  (J  -(^H  w  3  *  at  (2.16) 

and 

-U- 


’  %>V.V.VwV  -W. 


a  K  T 


(2.17) 


vm  = 

where  a  ts  a  pro por t i onal  i  t y  constant,  K  is  the 
Boltzmann  constant  and  T  is  the  electron  temperature. 

The  C^(y)  and  C^Cy)  are  the  Dingle  integrals  which  are 

7  7 

defined  as 


1  f"  xne_x 

Cn(y )  »  —  J  x-r— ~ T  dx  (2.18) 

o  J 

Using  these  new  collision  frequency  terms,  the 
susceptibility  tensor  M  takes  the  form  given  by  Eq. 
2.21,  which  is  equivalent  to  the  classical  form  of  the 
constitutive  relations  2.14  at  the  asymptotical  limits 
vm<<|  ui—  uihI  and  vm>>j  <*•* o>h|  with  ve  defined  as  the 
effective  collision  frequency  [32]. 


vm  <  < 


co 


ojh 


(2.19) 


\>m  2>  u  ♦  (uu 


(2.20) 


3 


2.3  -  WAVE  POLARIZATION  IN  THE  MAGNETOIONIC  MEDIUM 

The  polarization  of  a  plane  wave  propagating  in 
the  z-direotion  is  defined  as 


(2.22) 


In  a  magnetoionio  medium  the  following  relationships 
also  hold: 


(2-23’ 

If  B  0  is  chosen  to  lie  in  the  x- z  plane,  then  the 
constitutive  relations  2.13  become 


e0XEx=UPx+  i  nYPy 


e0  X  Ey  =  U  Py  -  inYPx+liYPz  (2.24) 


-  e0  X  Ez  =  U  Pz  -  i  HYP. 


Combining  Eqs .  2.22-2.24,2.25  leads  to  a  solution  of 
the  polarization  equation 


2  Y  ;  (  U  -  X  ) 


4  Y  2  L  ( ij  -  X  ) 


1  l  1  >2 


(2  OS'1 


where  YT  =  2,Y  and  Y^  =  nY  .  This  equation  shows  that 
there  exists  two  possible  pol ar iz ati ons  for  an 
electromagnetic  wave  propagating  in  a  given  direction 
in  a  magnet o i on i c  medium  [26]. 

2.4  -  THE  APPLETON-HARTREE  EQUATION 

Appleton  (1932)  and  Hartree  (1931)  developed  the 
theory  and  derived  the  refractive  index  for  vertically 
propagating  radiowaves  in  a  magnetoionic  medium  such  as 
the  ionosphere  [26,29,30].  Their  derivation  of  the 
refractive  index  is  often  referred  to  as  the 
Appi eton-Hartree  equation  which  is  given  by 


X 


U 


,.Y_T  + 

2  (  U  -  X  )  - 


CJ-X  )  2 


1  /  2 


(2.26) 


Eq .  2.26  demonstrates  the  complexity  of  radiowave 
propagation  in  the  presence  cf  a  superimposed  magnetic 
field.  Notice  that  in  general  there  are  two  possible 
refractive  indices  for  any  X,  due  to  the  ±  sign  in  the 
denominator.  This  shows  that  the  ionosphere  is  a 
doubly  refracting  medium.  Thus,  for  any  given 
direction  of  propagation  two  char  act er i s t i c  waves  are 
possible,  each  having  different  polarisations  and 
different  refract'  ve  indices.  These  two  characteristic 


waves  propagate  independently  and  have  different 
reflection  levels.  It  is  customary  to  call  these  two 
waves  "ordinary"  and  "extraordinary."  The  ordinary 
Co- )  wave  is  the  wave  component  whose  propagation 
character isti cs  are  influenced  minimally  by  the 
superimposed  magnetic  field.  On  the  other  hand,  the 
extraordinary  (x-)  wave's  propagation  is  influenced 
maximally  by  the  superimposed  magnetic  field. 

If  the  superimposed  magnetic  field  is  removed 
(Yl”Yt*°).  the  refractive  index  reduces  to  that  for  an 
isotropic  plasma. 


(2.27) 


If  collisions  are  neglected  and  the  wave  frequency, 
equals  the  plasma  frequency  (X  =  1),  n2  =  0.  This  is 
the  level  of  reflection.  Beyond  this  isvei  (X  >  1)  the 

wave  becomes  evanescent  and  does  not  propagate. 

The  A p pi  etc n-H ar t r ee  equation  exhibits  a  zero 
refractive  index  at  three  levels: 


X  -  1  ,  X  -  1  -  Y  ,  X  »  1  ♦  Y  (LT  -  1  ) 


(2.28) 


The  X 


1  reflection  level  for  the  o-wave  la  also  the 


reflection  level  in  an  isotropic  plasma.  X  »  1  -  Y  and 
X  »  1  +  Y  are  the  reflection  levels  for  x-waves.  There 

is  another  level  X  =  X^  where  n2  becomes  infinity  for 
a  collisionless  plasma,  or  very  large  when  the 
collision  frequency  is  small.  The  x-wave  may  be 
significantly  absorbed  at  this  level,  given  by 

Xr  -  f-r-ylz  (2.29) 

For  an  intermediate  inclination  of  the  magnetic 
field  and  Y  <  1,  the  variation  of  n2  with  X  is  as  shown 
in  Fig  2.1.  Figure  2.2  demonstrates  the  relative 
levels  of  reflection  and  absorption  in  terms  of  X  as 
well  as  the  regions  in  which  the  o-  and  x-waves  can 
propagate.  If  a  wave  is  incident  from  below  the 
ionosphere  with  a  monotoni ca 11 y  increasing  electron 
density,  the  x-mode  will  reflect  at  the  X  =  1  -  Y  level 
The  o-mode,  being  unaffected  at  this  level,  win 
propagate  up  to  X  -  1  and  reflect  there.  The  arrival 
times  of  these  two  reflected  components  will  be 
different  due  to  the  height  difference  between  these 
two  levels.  This  double  reflection  phenomenon  is 
called  "magneto- i oni c  splitting"  and  its  cause  is  the 
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Figure  2  2  Demonstration  of  mo d  e -  pr  o  p a ga  t  ion  regions  as 

a  function  of  the  plasma  frequency  parameter  X 


doubly  refracting  properties  of  the  plasma  containing 
the  Earth's  magnetic  field  [25,  26,  27], 

2.5  "  THE  Z-ECHO 

At  locations  where  the  Earth's  magnetic  field  is 
nearly  vertical,  experimenters  have  occasionally 
observed  a  third  echo  (named  z-echo)  during  ionosonde 
soundings.  Eckersley  (1950)  and  Rydbeck  (1950,1951) 
explained  this  observation  by  reflection  from  the  X  =  1+Y 
level  through  the  mechanism  of  mode  coupling  from  o-  to 
x-waves  at  the  X  =  1  level  [27].  If  the  angle  between 
the  wave  normal  and  the  Earth's  magnetic  field  is 
small,  n2  versus  X  behaves  as  shown  in  Fig.  2.3.  This 
figure  shows  that  around  X  =  1  the  curves  for  the 
o-mode  and  x-mcde  refractive  indices  are  nearly  equal. 
There  is  also  a  rapid  change  of  the  o-mode  polarization 
at  this  level.  This  rapid  change  of  polarization  may 
lead  to  a  coupling  between  the  o-  and  x-modes  and  cause 
the  transfer  of  energy  into  the  region  between  X  =  1 
and  X  »  1  +  Y  as  an  x-wave.  This  wave  would  then 
travel  up  to  the  level  X  =  1  *  Y  where  it  would  be 

reflected,  and,  after  coupling  back  to  o-mode 
polarisation,  it  would  appear  as  a  z-echo  on  an 
' one  gram  . 
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Ellis  pointed  oat  that  for  vertical  Incidence  the 
collision  frequency  in  the  F-region  of  the  ionosphere 
Is  not  great  enough  for  coupling  energy  effectively 
above  the  X  =  1  level  [35,36,37,38].  He  suggested  an 
alternative  explanation.  At  X  =  1,  If  the  wave  normal 
is  parallel  to  the  Earth's  magnetic  field  ( Y T  =  0), 
even  with  a  small  number  of  collisions  n  would  be 
non-zero,  which  may  cause  the  o-wave  to  penetrate  into 
that  region  through  mode  coupling. 

Fig.  2.4  shows  n2  versus  X  for  this  case.  It  is 
seen  that  the  refractive  indices  of  o-  and  x-modes  are 
almost  equal  at  X  =  1.  This  condition  may  always  be 
met  provided  the  o-wave  direction  of  propagation  is 
near  a  critical  angle  of  Incidence  0G  below  the 
ionosphere,  given  by 


sin  9, 


— — -  sin  9  a 

+  L 


(2.30) 


Having  penetrated  the  X  =  1  layer,  the  x-wave  then 
would  propagate  up  to  the  X  =■  1  +  Y.  level  and  would 

reflect  in  a  direction  which  would  not  be  observable 
from  the  ground.  However,  Ellis  surmised  that 
•'  r  r  e  gul  ar  i  t  i  es  at  the  X  =  1  +  Y  level  might  cause 

backsoattering  which  wend  then  be  observed  at  the 


transmitter  location  as  a  z-echo  [35].  It  should  be 
noted  that  the  term  z-echo  or  z-mode  wave  is  used  for 
historical  reasons.  Theory  shows  that  it  Is  actually  a 
downgoihg  .o-wave  which  underwent  successive 
t r ans f orm at i on  from  o-  to  x-  to  o-  polarizations  near 
the  X  -  1  level. 

2.6  -  THE  BOOKER  QUARTIC 

The  Appl  et  on-H  ar  t  r  ee  formulation  introduced  in 
the  previous  section  is  most  useful  for  vertically 
incident  waves.  For  oblique  incidence  a  more  general 
formulation  is  desirable  for  investigating  radio  wave 
propagation  in  a  magnetoionic  medium. 

Assume  a  wave  is  represented  by  the  following  wave 
function  in  which  the  y  dependence  is  not  included  for 
ease  of  discussion. 

exp(ik(ct  -  n(x  sin  9  +  z  cos  9))}  (2.30 

For  an  isotropic  medium  the  refractive  index  n  is 
independent  of  the  direction  of  the  wave  normal. 
Therefore,  to  determine  the  propagation  angle  9  of  the 
wave  for  a  particular  n,  Sneil’s  law  applies 


where  9 j  is  the  angle  of  incidence  in  free  space  (n=l) 
below  the  ionosphere. 


For  anisotropic  media  the  small  wave  function  can 
be  used  to  represent  either  o-  or  x-waves.  Sneii's  law 


still  applies,  giving 


sin  9j  =  nx  sin  6X  =■  nQ  sin  9C 


(2.33) 


where  (nx,9x)  and  ( nQ  ,  0O )  are  the  refractive  indices 
and  propagation  angles  of  the  x-and  o-waves  ,  ■  respec¬ 
tively.  However,  in  this  case,  because  of  the 
anisotropy  both  nx  and  n0  depend  on  the  direction  of 
the  wave  normal.  That  is,  nx  and  nQ  are  functions  of 
9X  and  90  .  Since  9's  and  n's  are  both  unknown,  one 
cannot  simpiy  use  Snell's  law  to  determine  the 
parameters  of  Eq.  2.31- 


m 


A  more  general  method  for  investigating  oblique 
propagation  was  developed  by  Booker  [38,39,26].  Booker 


f 


introduced  the  variable  q,  which  is  more  meaningful 


than  n  for  oblique  propagation  in  the  ionosphere. 
Consider  n  as  a  vector  in  the  direction  of  the  wave 
normal  inclined  at  an  angle  e  from  the  vertical,  as 
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1 


shown  in  Fig.  2.5.  Then  q  Is  defined  as  the  vertical 
component  of  n 


q 


Fig.  2.5  Relationships  between  the  propagation 
parameters  (n,  q,  6t, 

and  satisfies  the  following  relationships: 

sin  8j  =  n  sin  8  (  2  .  3  *0 

q  =  n  cos  9  (2.35) 

Substitution  of  these  variables  into  the  wave  function 
Eq.  2.31  gives 


e xp  j 1 k ( ct 


:<  sin  9  j  -  q  z  )  f 


(2.36) 


Booker  has  shown  that  q  is  a  root  of  a  fourth- 
order  polynomial  which  is  known  as  the  Booker  quartic: 

F(q)  =  aq1*  +  8q3  +  Yq2  +  Sq  +  e  =  0  (2.37) 

where 

a  =  U2(U-X)  -  U  Y  2  +  n 2  Y  2  X 
6  =  2nX  Y  (S^  =  S2m) 

Y  =  -2U(U-X)  (C2U  -  X)  +  2Y2(C2a  -  X)  (2.38) 

♦  XY  2  { 1 -C  2n  2+  (Sit  *  S  j  m  )  2  j 

<5  =  2C  2n  2XY  2  (S  i  t  *  S  2m  ) 

c  -  (U-X)  (C2U  -  X)2  -  C2Y2(C2U  -  X)  -  C2XY2(S,t*S2m)2 

The  four  roots  of  the  3coker  quartic  correspond  to 
the  upgoing  ordinary,  dcwngoing  ordinary,  upgoing 
extraordinary  and  downgoing  extraordi nary  waves.  The 
roots  are  in  general  complex  because  of  the  damping 
effect  of  collisions.  The  real  part  of  q  governs  the 
phase  change  of  the  wave,  and  the  imaginary  part 
governs  the  attenuation  of  the  wave  as  it  progresses  in 
the  medium.  Upgoing  and  downgoing  waves  can  be 
identified  from  the  imaginary  part  of  the  root.  The 
'conservation  cf  energy  principle  requires  that  the 


waves  attenuate  in  the  direction  of  energy  flow,  so  for 


a  wave  propagating  in  the  positive  z-direction 
(upgoing),  the  imaginary  part  of  q  must  be  negative. 
Conversely,  for  waves  travelling  in  the  negative 
z-direotion  (downgoing),  it  must  be  positive. 

In  certain  regions  of  the  ionosphere  it  is 
possible  that  the  roots  of  the  Booker  quartic  may  have 
real  and  imaginary  parts  with  different  signs.  This 
means  that  the  direction  of  energy  flow  and  phase 
propagation  are  different  [39,25].  This  is  explained 
by  the  fact  that  for  anisotropic  media  the  wave  normal 
and  ray  directions  are,  in  general,  different.  For  the 
same  reason  the  reflection  level  for  anisotropic  media 
is  the  level  where  two  roots  become  equal,  not  where 
q  =  0  . 

Fig.  2.6  shows  the  real  part  of  q  for  the  vertical 
incidence  case  when  c  her  e  are  no  collisions.  It  is 
seen  that  the  reflection  levels  for  the  o-  and  x-waves 
are  at  X  =  1  and  X  =  1  ±  Y,  respectively.  Also  note 

that  at  the  X  =  XR  level  the  x-wave  roots  approach 
infinity. 

2.7  -  RADIO  WINDOWS  AND  MODE  COUPLING  IN  THE  IONOSPHERE 
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Platteville  1=0 
Vertical  Incidence 


Fig.  2.6  Roots  of  the  Booker  quartic  for  the  vertical 

Incidence  case 


As  discussed  previously  in  Sect.  2.4,  an  o-w.ave 
may  penetrate  through  a  reflecting  layer  as  a  result  of 
mode  coupling  to  an  x-wave.  The  terms  radio  window  or 
holes  are  used  to  describe  the  occurrence  of  these 
possible  penetrations  [35.40,41], 

The  radio  window  of  interest  here  occurs  when  the 
propagation  is  in  the  magnetic  meridian  at  the  critical 
angle  of  incidence 

0G  =  ±  arcsin  ( l/ y— 7— 7)  (2.39) 

Figure  2.7  shows  the  sample  behavior  of  the  real  part 
of  the  roots  of  the  Booker  quartic  for  ec  <  0.  It  is 
seen  that  the  upgoing  branch  of  the  o-mode  and  x-mode 
is  continuous  at  the  X  =  1  level.  This  suggests 
coupling  from  o-  to  x-waves.  The  energy  coupled  to  the 
x-mode  will  propagate  up  to  the  X  =  1  +  Y  level  where 
it  i3  reflected  and  eventually  absorbed  near  the  X  R 
level  [25,28].  When  the  propagation  is  still  in  the 
magnetic  meridian  but  incident  at  an  angle  0C  >  0,  the 
real  part  of  the  roots  become  as  shown  in  Fig  2.8. 

Note  that  thi3  figure  is  symmetric  to  the  previous 
figure  about  the  x-axis.  In  this  case  the  upgoing 


c-wave  may  also  couple  to  an  upgoing  x-wave  at  the  X  *  1 


Roots  of  the  Booker  q 
incidence  (propagation 


level.  However,  this  x-wave  reflects  just  above  X  =  1 
(rather  than  at  X  =  1  +  Y)  and  is  eventually  absorted 
in  the  absorption  region  near  X  «  XR  [25,28]. 

Throughout  the  text,  the  Xr  level  is  assumed  to  be 
the  absorption  level  for  the  x-wave.  This  occurs 
because  of  the  fact  that  the  imaginary  part  of  q 
becomes  large  as  this  level  is  approached  when 
collisions  are  present.  Therefore,  in  this  region 
propagation  over  a  small  height  range  will  result  in 
absorption  of  large  amounts  of  energy  by  conversion  to 
heat.  For  high  power  radio  waves,  this  heating  may 
significantly  modify  the  plasma. 

An  important  related  equation  is  what  happens  to 
the  flow  of  energy  of  an  x-wave  near  the  Xr  level  when 
there  is  no  physical  mechanism  for  absorption?  This 
happens  when  there  are  no  collisions,  in  which  the 
non-e vanes  cent  imaginary  part  of  q  is  always  zero 
throughout  the  magnetoionic  medi  am.  It  has  been  found 
that  as  the  x-wave  approaches  the  Xr  level,  the 
direction  of  -energy  flow  becomes  horizontal.  This 
phenomenon  is  true  both  with  and  without  collisions, 
provided  the  collision  frequency  is  not  very  great 


Gi nzbu 


The  x-z  plane  is  chosen'  to  be  the  plane  of  incidence, 
and  plasma  stratification  is  assumed  to  be  in  the 
z-direction.  The  incident  plane  wave  is  taken  to  be 
inclined  at  an  angle  0  from  the  z-axis.  The 
inhomogenous  medium  may  be  thought  of  as  being  made  ip 
of  a  very  large  number  of  discrete  strata  in  each  of 
which  the  electron  density  and  the  collision  frequency 
are  constant,  that  is,  homogenous.  Then,  in  each 
strata  the  field  quantities  vary  with  x  as  exp(ikxsin 
9[  and  do  not  vary  with  y,  thus  the  differential 
o  per  ators 

3  ,  3  ,  3  are  equivalent  to  -iksine,  0,  d/dz.  If  the 

3X  3  Y  3  Z 

factor  expjikxsin  0}  is  dropped  as  was  done  with 
exp{lwt},  the  equations  are  functions  of  z  only  and  the 
problem  becomes  one  dimensional. 

Combining  Maxwell's  curl  equations  with  the 
constitutive  relations  and  eliminating  Ez  and  Hz 
results  in  four  simultaneous  linear  differential 
equations  of  first  order  with  Ex,  Ey,  Hx,  Hy  as 
dependent  variables  [45,32]. 


u  iiTLW  AAAn  w 


■—  Ex 

dz  x 


dz  Ey 


ik|-Hy  +  ^Tm  (sineHy  +  Mzx^x 


*  ik  H, 


M  z  y  ^  y  )  1 


—  Hx  =  lk(MyXEx  +  (I  +  Myy)Ey  -  j  +  (sin8Hy  MZXEX  + 


MZyEy)  -  5  i  n  2  9  Ey 


(  2  .  H  1  ) 


—  Hy  =  i  k  { -  ( I  MXX)EX  -  HXyEy  -  (sin@H, 


^zx^x  +  Mzy 1 


Ez  and  Hz  are  given  as 


(sin0  y  +  MzxEx  +  MzyEy} 


( 2. U2) 


H7  =  sine  E , 


In  general,  in  an  homogenous  medi  im  these  four 
simultaneous  linear  differential  equations  result  in 
four  independent  solutions  for  the  field  quantities  Ex, 
Ey,  Hx ,  Hy .  These  four  independent  solutions 
correspond  to  the  upgoing  ordinary,  downgoing  ordinary, 
upgoing  extraordinary  and  downgoing  extraordinary 
waves,  propagating  independently.  Under  most 
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circumstances  ,  separation  of  the  wave  fields  in  this 
manner  is  also  applicable  in  inhomogenous  media. 
Exceptions  to  this  occur  at  the  levels  of  reflection 
and/or  coupling  where  the  independence  of  the  solutions 
are  lost. 

Eqs  .  2.41  may  be  written  in  matrix  form  as: 


~  (a)  -  ikT  (e) 


where  e  is  a  4  x  1  column  vector 


e  *  (Ex  Ey  Hx  Hy)T 


and  T  is  a  4  x  4  matrix  given  by 


a i n0M, 


3lnlM7, 

1  ♦  M, 


’y*  i  ♦  M. 


7.'  7.X  u  ,n 

1  »  M , _  ■  'XX  1  ' 

z  z 


cos9  ♦  M 


yy  '  ,  m 


»X7»7.i 

1  ♦  M , 
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(2.43) 


(2.4  4)- 


0  - (  003  2  ,  *  M, T  ) 

1  *  M , , 


(2.45) 


■  3  i  n  9  M ,, 

T-.  Tit 


3  1  n , Mv ■ 

1  *•  M, . 


In  each  stratum  the  field  components  depend  on  z 

only  through  the  factor  exp[-ikqz[  which  gives 

d  5  -ikq.  Eq.  2.43  then  can  be  written  as 
d  z 

Te  =  -qe  (2.46) 

The  nontrivial  solutions  for  q  are  the  eigenvalues  of 
the  T  matrix  or,  equivalently,  the  solution  to  Eq.  2.47 

det (T+ql  )  =  0  (2.47) 

which  is  another  form  of  the  Booker  quartic.  Expan  ' ng 
the  determinant  gives  the  coefficients  of  the  quartic 

o  =  1 

S  *  Tii  *  T4  4 

y  =  T  1  1  T  4  4  -  T  4  T  u  1  *  T^2 

5  =  -T32(T,  ,  +  T  ^  4  )  +  T  3  4  T  4  2  +  T i  2  T  3 1 
e  =  ^  1  1  (  T  3  4  T  4  2  “  T44T32)  +  1  2  ^  3  1  T  4  4  _  T  3  4  T  4  ) 

+  T14^T32T41  "  T  3 i T  4  2  ^ 

For  each  eigenvalue  there  is  an  eigenvector  e  1 
which  satisfies  the  equation 


solutions  are  not  considered.  Rather,  the  exact 


solution  of  the  field  components  is  required 
necessitating  use  of  the  full  wave  solution  technique. 

The  full  wave  problem  is  to  find  the  exact 
solution  to  the  differential  Eqs  .  2.41  either 
analytically  or  numerically.  Analytical  solutions  are 
limited  to  special  kinds  of  electron  density  profiles 
and  magnetic  field  conditions  [26].  Numerical 
solutions,  on  the  other  hand,  are  more  general  and  are 
therefore  used  here  [32,44]. 

The  differential  equations  2.41  have  four  sets  cf 
Independent  solutions  and  require  four  sets  of  boundary 
conditions  for  their  complete  solution.  At  the  bottom 
level  where  the  transmitter  is  located  the  upgoing 
incident  wave,  which  gives  two  of  the  boundary 
conditions,  is  assumed  to  be  known.  However,  the  other 
two  resulting  from  dcwngcing  waves  are  net  known,  which 
implies  that  boundary  conditions  at  the  bottom  x  e  ve  a 
are  insufficient.  Above  the  highest  reflection  level 
only  upgoing  waves,  which  supply  the  other  two  boundary 
conditions,  can  exist.  Here  upgoing  waves  are  defined 
to  be  waves  travelling  in  the  positive  z-direction  ,  or 
are  evanescent  with  decreasing  amplitudes  in  the 
positive  z-direction.  The  down  go  i  ng  wave  com  po  n  en.es 
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above  the  highest  reflection  level  mast  be  zero;  this 
completes  the  required  four  boundary  conditions.  Two 
sets  of  solutions  are  obtained  from  these  upper 
boundary  conditions  by  separately  integrating  them  down 
through  the  ionospheric  model.  Since  these  two 
solutions  are  independent,  the  total  wave  field  is  then 
a  linear  combination  of  them.  The  coefficients 
relating  the  two  solutions  to  the  total  wave  field  is 
found  by  extracting  their  upgoing  components  and 
matching  them  to  the  known  incident  wave  at  the  bottom 
level  where  the  integration  is  ended. 


2.10  -  DERIVING  THE  BOUNDARY  CONDITIONS 

At  the  highest  level  where  integration  is 
started  the  initial  values  of  the  upgoing  field 
components  Ex,  Ey,  Hx,  Hy  are  required  as  boundary 
conditions.  The  mediim  above  the  highest  level  is 
considered  homogenous  so  that  the  field  quantities  vary 
with  z  as  -ikq  as  before.  Substituting  this  into  the 
Eq3 .  2.41  gives  the  wave  polarizations  for  each  field 
component  [32] 


-  M  Y  •/  (  co  s  2  9  ♦  M ..  )  +  M  (  q  s  i  n  0  +  M  v  ) 


'  1  -  q  *  xx  )  (  cos  '  9  *•  M  2  z  ;  -  (  q  s  i.  n  3  M  2  x  )  (  qsin-3  +  M 
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- 1 


I 


zz 


( si n6Hy>M2xnx  +  MZy  ) 


(2.52) 


H 


2  =  sineEy 


The  two  boundary  conditions  correspond  to  the 
quartic  roots  with  negative  imaginary  parts  which 
denote  upgct.ng  waves.  Solving  the  quartic  and 
substituting  appropriate  roots  into  Eq.  2.51  results  in 
two  sets  of  initial  solutions. 


1 1 


=  (  E  1  v  E  1 


X  -  y  H'X  H‘y)T 


(  E  2  y  E2  H  2  y  H2V)T 


(2.53) 


(2.5*0 


2.11  -  INTEGRATION 

Numerical  integration  is  done  by  the  fourth- 
order  Runge-Kutta  method.  Integration  starts  with 
initial  solutions  (1)  and  (2)  at  the  maximum  height  and 
proceeds  downwards  in  the  medium  with  h  meter 
decrements.  At  selected  intervals  cf  h  ,  s  c  i  it  >  c  ns  '  1  ) 
and  2  )  are  tore  d ;  t  n  e  v  w  >  „  i  - 


ise  d 


U  ar 
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reconstruct  the  total  field  components  by  forming  their 
appropriate  linear  combination.  The  Runge-Kutta  scheme 
is  outl i ned  below  . 

The  simultaneous  differential  equations  can  be 
represented  as 


y'i  -  f i  (y i ,  y 2 .  y  3 .  y*.  z ) 


where 


1=1,  ,  3,  4 

(2.55) 


Y\  =  Ex  y2  =  Ey  y3  ~  Hx  y*  =  Hy  (2.56) 


Then  the  integration  coefficients  become 


+  1  /2k 


*1  /2k 


kii  =  -h  fj.  (yi  ■  y 2 '  ""h.  y» 

v  +  1  /2k  7-U2h) 

J  4  4  0  »  **  / 


+  1  /2k 


+  1  /2k 


+  1/2k  v  +  1  /  2  k 


ki  2  =  -h  f  i  (y,  ''“n,  y2  2  i  *  :  3  3i* 

y/172kM,  z~  1  /2h  ) 

kl  ,  =  -h  f •  (v  +1/2k  v  +  1/2k  v  +  1/2k 

Ki  3  12*  J  2  2  2  *  )  3  32* 

v  + 1 /2k  7- 1 /2h  \ 

y  4  4  2  »  «  1 

k  1 4  =  -  h  (y,+  k13,  y-2+  k23,  y}+  k33,  y4+  k,3,  z  -  h  ) 
y  i  ( z-  h )  =*  y  t  ( z )  +  ^  (kla  +  2\<il  +  2klz  +  ki3) 


2.12  -  NUMERICAL  SWAMPING 

It  has  been  observed  that  daring  integration 
numerical  swamping  might  cause  problems  because  of  the 
limited  accuracy  of  the  computations.  For  example,  at 
the  beginning  of  integration  when  the  initial  solutions 
are  calculated,  or  during  integration,  numerical  errcrs 
such  as  roundoff  may  introduce  a  small  fraction  of  an 
evanescent  wave  into  the  travelling  wave  solution. 

Thi3  evanescent  wave,  which  is  very  insignificant  at 
the  beginning,  may  eventually  dominate  the  travelling 
solution  and  lead  to  incorrect  results. 

Several  methods  have  been  devised  to  overcome  this 


pro  bl em  . 


One  method  is  tc  change  dependent  variables 


and  transform  the  differential  Eq.  2.43  to  a  different 
form.  Instead  of  Ex,  Ey,  x,  y,  four  elements  of  the 
admittance  matrix  A  may  be  used  [40] 


Alternatively,  elements  of  the  variable  p  are  also 
poss i bi e  [ 40  ] 


p  -  (A  -  I)  (A  +  I ) ~ 1  (2.59) 

Since  solutions  (1)  and  (21)  are  both  contained  in 
the  wave  admittance  matrix  A,  these  two  methods  require 
only  one  integration.  However,  they  yield  only  the 
ratios  of  field  components  rather  than  the  true  field 
amplitudes,  and  are  employed  primarily  for  calculating 
reflection  coefficients.  If  the  field  variation  with 
height  is  required,  the  choice  of  independent  variables 
must  in  each  case  yield  Ex,  Ey,  Hx,  Hy  and  the 
differential  Eq.  2.43  must  retain  their  original  form. 

The  namericai  swamping  problem  has  been  overcome 


by  a  different  method  sjggested  by  Pitteway  [ n  4  ]  . 


T  he 


method  is  based  on  constraining  solutions  at  regular 


intervals  during  integration  and  forcing  them  to  remain 
independent.  This  is  accomplished  by  adding  a  fraction 
of  one  solution  to  the  other  one.  This  is  possible 
because  the  differential  equations  are  linear.  In 
practice,  these  modifications  are  recorded  along  with 
the  field  components  during  integration  and  can  be 
accounted  for  later,  during  the  reconstr uction  of  the 
total  field. 

The  fraction  of  one  solution  added  to  the  other 
solution  must  satisfy  certain  conditions.  Let  e2 
represent  the  solution  to  be  constrained 


(  E  2  x  E2y  H 2  x  H 2  y  )  T 


and  e'2  the  solution  to  replace  e 


e '  2  =  e  2  +  ae, 


(2.60) 


(2.61  ) 


where  a  is  a  complex  constant  such  that  e'2  is 
Hermitian  orthogonal  to  e,.  That  is,  the  dot  product 


e '  2  •  e  *  ,  =  0 


(2.62) 
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where  *  denotes  the  complex  conjugate.  This  requires 
that 


e*  i  *  !i 


(2.63) 


2.13  "  MATCHIM  G  SOLUTIONS  TO  THE  INCIDENT  WAVE 

Below  the  Ionosphere  where  the  medium  Is  assumed 
to  be  homogenous,  the  known  incident  wave  must  be 
matched  to  the  upgoing  wave  components  of  both 
solutions  through  an  appropriate  linear  combination. 

If  the  incident  wave  components  are  Exinc,  Eyinc  and 
the  upgoing  field  solutions  are  Ex‘,  Ey  1  ,  Ex2,  Ey2, 
then 


EYinc  =  aE  1  x  +  S  E  2 ■ 


Eyinc  =  a  E 1 y  ♦  S  E  2  y 


(2.64) 


[  ->  so 


Solving  for  ct  and  S  gives 


F  i nfir  2  -  r  1 ncp  2 

Dx  cy  Ly 

f  r  f  Z  .  "p  1  p  2 
HXC,y  C,yCX 


(2.66) 
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EyincElx  -  ExincE‘y 
E  1  x  ^  2  y-  ”  E  1  2  x 


(2.67) 


Another  method  of  mat  chi  ng-  e  xi  sts  when  the 
incident  wave  has  only  o-  or  x-wave  polarization.  Let 
£y,0inc>  £y,x*nc  denote  the  o-  and  x-wave  components, 
respectively.  Then, 


Ey,Q^nC  =  ClE  ‘  y  (  Q-  +  8E2yj0 


Ey  ,  X  1  nC  =  aE  1  y  (  x  +  B  E  2  y  f  X 


(2.68) 


(2.69) 


If  there  is  only  an  apgoing  o-wave,  then  Ey>xinc  =  0 
and  the  coefficients  a  and  8  become 


a  =  -8 


c  i  n  c  tr  i 
ty  ,o  L  y  ,  x 


3  =  pi 

t  y , 0L  y , x 


^  y , O1  y  ,  x 


(2.70) 


(2.71) 


The  coefficients  for  the  x-wave  can  be  derived  in  a 
s i mi  1  ar  manner  . 

After  the  coefficients  a  and  8  are  determined,  the 
total  fields  anywhere  in  the  ionosphere  can  be  found  by 
taking  the  linear  combination  of  the  two  solutions 
which  were  stored  d  jr ' ng  integration. 
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2.14  -  DECOMPOSING  TOTAL  FIELD  COMPONENTS  IN  THE 
LOWER  BOUNDARY 

The  fall  wave  solution  technique  gives  the  total 
fields.  One  would  also  want  to  know  the  upgoing  and 
downgoing  o-  and  x-wave  fields  for  determining  the 
reflection  properties  of  the  ionospheric  model. 

Decom posing  the  total  fields  into  their  magnetoicnic 
modes  in  an  homogenous  medium  is  possible  using  the 
method  described  below. 

Let  subscripts  i  =  1,2,3,4  refer  to  the  upgoing  c-, 
downgoing  o- ,  upgoing  x-  ,  downgoing  x-waves, 
respectively.  The  total  field  may  then  be  written  as 
a  sum  of  these  four  waves  such  that 

E  I  St  (2.72) 

H  l  H*  (2.73) 

where  only  the  total  fields  E  and  H  are  known.  These 
equations  may  be  written  in  terms  of  the  polarization 
states  o  of  each  wave  component  (Eqs.  2.75-2.76)  where 


(2.74) 


2.51-2.52  with  its  respective  root  q.  Now,  the  total 
fields  (Ex  Ey  Hx  Hy)T  on  the  left  of  Eq.  2.75  and  the 
4x4  polarization  matrix  are  known.  Solving  Eq.  2.75 
for  (Ey,  E y 2  Ey,  Ey„)T  gives  the  jpgoing  o- ,  downgoing 
o- ,  jpgoing  x- ,  downgoing  x-wave  components  of  Ey, 
respectively.  The  other  cartesian  components  may  then 


be  fcjnd  from 
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CHAPTER  3 

MODEL  COMPUTATIONS 

In  this  chapter  field  solutions  for  simple 
ionospheric  models  are  performed  using  the  full  wave 
solution  program  developed  by  Seliga  [32].  The 
solutions  give  the  wave  fields  for  an  incident  o-wave 
for  different  electron  density  profiles,  collision 
frequencies  and  angles  of  incidence.  With  the  aid  of 
the  corresponding  Booker  quartic  solutions,  the 
behavior  of  electromagnetic  waves  in  the  plasma  are 
determined.  The  reflection  coefficients  for  several 
propagation  directions  are  also  calculated  to  determine 
the  shape  and  angular  width  of  the  radio  windows.  The 
results  demonstrate  the  importance  of  o-  to  x-mode  wave 
coupling  which  occurs  when  propagation  is  near  the 
critical  angle  in  the  magnetic  meridian. 

3.1  -  PARAMETERS  OF  THE  IONOSPHERIC  MODEL 

Three  different  models  are  chosen  which 
approximate  the  increasing  electron  density  of  the  F 


region  of  the  ionosphere. 


Electron  Density  Profile 


Each  profile  has  a  linearly  increasing  electron 
density  with  height  given  by  the  equation 


N(z)  =  NB  *  NS  (-) 


where  Ng(m  3)  is  the  initial  electron  density  at  the 
reference  height 


z  =>  0 

Ng/H  (m-1*)  is  the  rate  of  change  of  electrons 
density  with  height 

z  (m)  is  the  height  relative  to  the  reference 


level . 

The  three  models  chosen  are  shown  in  Table  3.1 


Table  3 • 1  -  Parameters  of  the  selected  ionospheric  models 


Model  1 


Model  2  &  3 


Ng  ( m- 3 ) 


1.8  x  10 


1.8  x  1 0 1 1 


Ns  ( m“ 3 ) 


2.7  x  10 


2.7  x  10 


H  (m) 
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Frequency  : 


All  computations  are  performed  at  a  wave  frequency 
5  MHz 


Collision  Frequency  : 

A  constant  collision  rrequency  is  used  throughout 
the  medium.  For  each  model,  vm  <<  |  u  -  |  so  that  the 

effective  collision  frequency  ve  is  given  by  Eq.  2.19. 
Table  3.2  gives  the  collision  frequency  values  for  each 
model. 


Table  3 • 2  -  Collision  frequencies  for  the  ionospheric  models 


Model  1 


Model  2 


Model  3 


vm  (3  1  ) 


ve  (  s"  1  )  2.5  x  10 


Geomagnetic  Field : 


2,5  x  10 


The  geomagnetic  field  conditions  at  Platteville, 
Colorado  (Fig.  3-1)  are  chosen,  since  a  major  heating 
facility  was  located  there  and  much  experimental  data 
exists  at  this  location. 
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Table  3.3  -  Critical  level  heights  for  the  ionospheric  models 

Model  1  Models  2  &  3 

X  =  1  48 1 8  m . 

X  =  1  ♦  ¥  8447  m. 

X  =  1  -  Y  11 90  m. 

X  =  XR  4639  m. 

Critical  Angle  of  Incidence : 

The  critical  coupling  angles  for  each  model  are  same 

9  q  =  ±  10.696  ° 

3.2  -  PARAMETERS  OF  THE  FULL  WAVE  SOLUTION  PROGRAM 

Cartesian  components  of  the  E  and  H  fields  are 
outputted  by  the  full  wave  solution  program.  Field 
magnitudes  are  relative  to  the  incident  wave  amplitude. 

Note  that  since  propagation  is  in  the  x-z  plane,  Ey  and 
Hy  are  the  components  per pendi c ul ar  to  the  plane  of 
propagation. 

Integration  Increment  : 


481.8  m  . 
84  477  frT. 
119.0  m . 
463  .  9  m  . 


The  integration  increment  h  is  determined  by  a 
computational  procedure.  Three  arbitrary  angles  are 
chosen  for  each  model,  one  in  the  S-N  direction  near 
the  critical  coupling  angle,  one  in  the  vertical 
direction  and  one  in  the  N-S  direction  near  the 
critical  coupling  angle.  For  each  of  these  angles  the 
program  is  run  several  times,  each  with  decreasing 
integration  increments  until  the  reflection 
coefficients  and  fields  converge.  Then  the  program  is 
run  for  the  rest  of  the  angles  using  the  resulting 
integration  increment.  Table  3.^  shows  the  integration 
increments  used  for  the  models. 


Table  3.4  -  Integration  increments  used  for  each 

ionospheric  model 


Model 

3 

1  0  3 

m  . 

200 

s  "  1 

0 

1  0  3 

m  . 

0 . 02 

m  . 

Model 

2 

1  0  3 

m  . 

1  0  3 

3“  1 

0 

1  0  3 

m  . 

0  .  1 

m  . 

Model 

1 

(N-S) 

1  0  * 

m . 

1  0  3 

3  ~  1 

1  0  3 

6x  1  0  3 

m  . 

0 . 05 

m  . 

Model 

1 

(  S-N) 

1  0  “ 

m . 

1  0  3 

3“  1 

1  0  3 

9x  1  0  3 

m . 

0.05 

m  . 

Hermitian  Orthogonal  Constraint : 

-  C-osve-fe-r a i n i n g  of  the  solutions  is  done  at  every  2  0  to 

50  integration  steps,  depending  on  the  storage 
requirements  of  the  computer  program. 
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Computation  Time: 

The  program  was  run  on  an  IBM  3081  with  the  VS 
FORTRAN  compiler.  The  average  run  times  for  the  models 
are  given  in  Table  3.5 

Table  3 . 5  -  Average  execution  times  for  the 
ionospheric  models 

Total 

Execution  Time  Number  of  An gles  Time 


Mode  1 

1 

(  S-N) 

9  min/angle 

22 

1  98 

min. 

Model 

1 

(  N-S) 

6  min/angle 

22 

1  32 

min. 

Model 

2 

1  min/angle 

60 

60 

min. 

Model 

3 

3  in/angle 

60 

1  80 

min. 

3-3  -  INTERPRETATION  OF  THE  RESULTS 

In  this  section  the  field  solutions  for  Model  1 
at  9  -  0  and  9  =*  ±  8C  in  the  magnetic  meridian  are 
presented  along  with  the  corresponding  roots  of  the 
Booker  quartic.  The  levels  X  »  1 ,  M  I,  1  -  K,  Xp  are 
indicated  on  the  full  wave  solutions  as  vertical  lines 
numbered  1,  2 ,  3,  4 ,  respectively.  On  top  of  tne 
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A  /.  ,\V, 


figures  are  the  parameters  of  the  ionospheric  model, 


including  Ng,  Ng,  H,  vm  which  are  denoted  by  C0N(1), 
C0N(2),  C0N(4),  C0N(5),  respectively.  The  solutions 
for  other  models  and  magnetic  fields  are  given  in 
Appendix  A. 


i 


Case  1  :  Model  1,  Vertical  incidence  (8  =  0) 

Figures  3.2,  3.3  and  3.4  show  the  Ex,  Ey ,  E2 
components  of  the  vertical  incidence  solution, 
respectively.  The  upgoing  incident  wave  has  ordinary 
polarization.  The  standing  waves  below  X  =  1  are  due 
to  interference  between  the  incident  upgoing  o-wave  and 
reflected  o-wave.  The  reflection  coefficients  for  E 
and  H  are 


i 


R( E)  =  0.889 
R(  H)  =  0  .  389 


I 


which  indicate  that  most  of  the  incident  energy  is 
reflected  back  at  vertical  incidence. 

The  roots  of  the  Booker  quartic  are  useful  for 
understanding  the  observed  reflection  mechanism.  Figs. 
3.5-3. 5  show  the  real  and  imaginary  parts  of  the  roots. 
The  roots  corresponding  to  the  upgoing  o-wave  and 
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downgoing  o-wave  are  denoted  by  A  and  B,  respectively. 
The  reflection  level  for  the  o-wave  occurs  where  the 
two  roots  become  equal,  which  is  X  =  1  in  this  case. 

The  imaginary  part  of  the  roots  cor r es pond i ng  to  A  and 
B  appears  to  be  zero,  indicating  that  there  is  no 
energy  lost  during  this  process.  Actually,  q  has  a 
small  imaginary  part.  A  detailed  plot  would  reveal 
that  Im(q)  for  A  is  negative,  indicating  that  direction 
of  energy  flow  is  upwards;  Im(q)  for  B  is  positive, 
indicating  that  direction  of  energy  flow  is  downwards. 

A  sketch  of  the  ray  path  for  the  vertical  incidence 
case  is  shown  in  Fig.  3-7.  The  rays  bend  towards  the 
north  as  the  X  »  1  level  is  reached,  and  become 
per  pend i cul ar  to  BQ  at  that  level.  More  detailed 
information  on  these  ray  paths  may  be  found  in 
[26 ,42,43]. 

Beyond  X  =  1,  Re(Q)  =  0  and  Im(q)  <  0  for  the  o-wave 
as  denoted  by  C  in  Fig.  3.6.  This  indicates  that  the 
o-wave  is  evanescent  and  does  not  propagate  in  this 
region.  In  the  full  wave  solutions  this  effect  is  seen 
as  a  sharp  decrease  of  the  field  amplitudes  beyond  X  = 


Another  interesting  property  is  the  existence  of  an 
Sz  component,  although  the  propagation  vector  i  a  in  th- 
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z-direction.  This  is  common  to  anisotropic  media.  The 
magnitude  of  Ez  is  quite  insignificant  at  the  bottom 
level  (Fig.  3.4),  but  as  the  waves  progress  up  to  the 
X  =»  1  level  it  grows  larger. 

Another  observation  of  interest  is  th-e  changing 
wavelength  of  the  fields,  decreasing  with  height.  This 
is  in  agreement  with  the  roots  of  the  quartic  (equal  to 
the  refractive  index  in  this  case)  which  indicate  that 
n  ♦  0  as  X  •*  1  . 

Case  2 :  Model  1,  Critical  Angle  of  incidence  (S-N): 

Propagation  Towards  Magnetic  North 

Figures  3.8,  3.9  and  3.10  show  the  Ex,  Ey ,  Ez 
components  of  the  fields,  respectively.  The  field 
magnitudes  in  this  case  are  nearly  constant  at  the 
bottom  level.  The  absence  of  any  standing  wave  peaks 
below  X  R  indicate  the  absence  of  significant  reflected 
downgoing  waves.  The  reflection  coefficients  for  E  and 
H  are 


R( E)  =  6 . 4450  x  1 0" 9 


'TV  « J  .JJ  KTTTOI 


R(  H) 


6.9295  x  10“3 


Wave  Frequency  :  5.0  E  6  Angle  of  Incidence 

Magnetic  Y  Vector  :  3.150  E-1  Con  (1) 

Tot.  Geomag.  Field  :  5 . 643  E-5  Con  (2) 


which  implies  that  practically  all  of  the  incident  wave 
is  absorbed  in  the  medium. 

Figs.  3-11  and  3-12  show  the  real  and  imaginary 
parts  of  the  Booker  quartic  roots.  The  upgoing  o-wave 
path  is  denoted  by  A.  Unlike  the  vertical  incidence 
case,  the  roots  do  not  follow  path  B  corresponding  to 
the  downgoing  o-wave,  but  progress  up  to  just  above 
X  *  1  where  the  two  roots  are  equal,  then  become  very 
large  at  the  XR  level.  This  is  more  clearly  seen  in 
Figs.  3-13  and  3 • 1 4 . 

The  roots  of  the  3ooker  quartic  show  that  upgoing 
o-waves  (branch  A)  propagate  up  to  just  above  X  =  1 
(X  *  1.0015  for  this  case),  reflect  back  and  approach 
X  R  .  The  corresponding  ray  path  for  this  case  is 
sketched  in  Fig.  3*15.  The  phenomenon  described  above 
is  seen  in  the  field  solutions  as  standing  waves  above 
the  X  R  level.  Above  X  =  1.0015  the  waves  exnibit 
evanescent  behavior  and  decay  rapidly  as  expected  from 
the  quartic  roots. 

The  imaginary  part  of  q  (Fig.  3.14)  shows  that  Im(q) 
is  zero  or  almost  zero  for  the  upgoing  o-wave  (branch 
A).  On  the  way  back,  after  reflection  as  an  x-wave, 

I M  (  q  ;  becomes  large  near  X  a .  This  increase  becomes 
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seen  by  the  following  branch  C  in  Fig.  3  •  1  **  •  The 
results  show  that  for  this  critical  incidence  case 
(S-N),  most  of  the  incident  energy  is  focused  and 
absorbed  in  a  very  narrow  region  just  below  the  X  =  1 
level. 

Near  the  Xp  level  the  real  part  of  q  along  C  becomes 
very  large  and  Ez  exhibits  a  peak  of  40dB  compared  to 
an  incident  wave  intensity  of-  around  5dB.  This 
indicates  that  E  and  k  become  almost  vertical  at  Xp 
while  the  direction  of  energy  flow  is  horizontal. 


Case  3 :  Model  1,  Critical  Angle  of  incidence  (N-S): 

Propagation  Towards  Magnetic  South 
Figures  3-16,  3*17,  and  3*18  show  the  Ex ,  Ey ,  Ez 
components  of  the  fields,  respectively.  The  field 
magnitudes  in  this  case  are  also  nearly  constant  at  the 
bottom  level.  The  absence  of  standing  waves  below  X  p 
indicates  the  absence  of  significant  reflected 


downgoing  waves.  The  reflection  coefficients  for  E  and 
H  are 


R  (  E )  =  6.3  x  1  0"  3 


R ( H )  =  o.3  x  10' 
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which  indicate  that  almost  100$  of  the  incident  energy 
is  absorbed  in  the  medium. 


The  roots  of  the  Booker  quartic  (Figs.  3.19-3-20), 
corresponding  to  the  upgoing  o-wave,  is  continuous  at 
the  X  =  1  level.  The  roots  continue  (branch  B)  up  to 
the  1  +  Y  level,  return  back,  with  the  real  part 
getting  large  as  the  Xr  level  is  approached  (branch  C). 
These  indicate  that  the  upgoing  o-wave  progresses  up  to 
the  X  »  I  level,  couples  to  the  upgoing  x-wave  and  is 
reflected  back  from  the  1  +  Y  level.  This  is  seen  on 

the  field  solutions  as  standing  waves  between  XR  and 
1  +  Y.  The  reflected  x-waves  do  not  attenuate 

significantly  until  the  Xr  level  is  reached;  this  is 
observed  on  the  roots  of  the  Booker  quartic  as 
imaginary  parts  being  zero  (Figs.  3.20,  3.22,  branches 
A ,  B,  C).  Figures  3-20  and  3-21  show  that  the  roots 
cross  the  Xr  level  (branch  D),  contrary  to  the  3-N 
propagation  case  in  which  the  roots  reach  their  maximum 
at  the  Xr  level . 

The  imaginary  part  of  q  (Fig.  3.22,  branch  D)  is 
large  between  Xr  and  1  -  Y.  The  waves  attenuate 
significantly  and  the  energy  is  absorbed  in  this 
region,  most  of  it  localized  near  XR.  Fig.  3.23  shows 
the  s/. etch  of  tne  ray  path  for  this  case.  Note  that 
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near  Xr,  the  real  part  of  q  and  Ez  component  is  very 
large.  This  indicates  that  E  and  k  are  almost'  vertical 
while  the  energy  flow  is  horizontal. 

3.4  -  WIDTH  OF  THE  COUPLING  WINDOW 

There  is  a  region  around  the  critical  coupling 
angle  at  which  the  o-wave  to  x-wave  coupling  is  most 
effective.  The  width  of  this  region,  which  is  referred 
to  as  the  coupling  window,  is  determined  by  performing 
a  series  of  full  wave  solutions  near  the  critical  angle 
of  incidence.  The  incidence  angle  is  varied  around  the 
coupling  angle  in  two  planes:  (i)  in  the  magnetic 
meridian  (N-s),  and  (ii)  perpendicular  to  the  magnetic 
meridian  (E-W).  The  -  3  d  B  points  are  used  to  define  the 
edges  of  the  coupling  window.  Table  3.6  lists  the 
angular  widths  for  the  E  field.  Figures  3.24-3.29  and 
3-30-3.35  show  the  variation  of  the  reflection 
coefficients  with  the  angle  of  incidence,  in  the 
magnetic  meridian  and  off  the  magnetic  meridian, 
resnecti vely .  The  rest  of  the  reflection  coefficients 
are  listed  in  Appendix  B. 

The  results  show  that  the  angular  shape  of  the 
coupling  winders  in  a  0  cro  x  i  ma  t  e  1  v  circular  around  the¬ 
oretical  coupling  angle.  The  gradient  a  f  t  n  e  electron 
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Table  3 • 6  -  Angular  Widths  of  the  Coupling  Windows 


Model 


N-S  and  S-N 


Angle  -3  dB ( low )  ~3dB ( h  igh  )  Width  -3  dB 
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FLg.  3-27  Reflection  coefficient  vs.  angle  of 
incidence:  Model  2,  in  the  magnetic  meridian,  toward 

magnetic  3  o  u  t  h 
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density  appears  to  be  a  major  determining  factor 
affecting  the  width  of  the  coupling  window.  The 
angular  width  of  Model  1  is  about  3-5-^°  and  Model  2  is 
about  10-11°.  The  parameters  of  Model  1  and  Model  2 
differ  only  in  the  gradient  of  the  electron  density  H 
(H  Model  1  H  =  1 0 - ,  Model  2  H  >  10*).  This  indicates 
that  steep  gradients  are  more  favorable  for  effective 
coupling.  Also,  for  a  given  ionospheric  model,  the 
coupling  window  in  the  S-N  direction  (+10. 7°)  appears 
to  be  wider  than  the  window  in  the  N-S  direction 
(-10.7°)  by  about  0.5  -  1 0  . 

No  significant  difference  is  found  between  Model  2 
and  Model  3  which  differ  only  by  the  collision 
frequency.  The  collision  frequencies  in  these  models 
do  not  seem  to  be  a  determining  factor  affecting  the 
width  of  the  coupling  window,  or  are  of  secondary 
importance  compared  to  the  gradient  of  electron 


density. 


CHAPTER  4 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  this  research  was  to  demonstrate 
the  possible  role  of  o-  to  x-  mode  coupling  in 
ionospheric  radiowave  heating  experiments.  Most, 
previous  related  studies  and  experiments  have  not 
considered  the  importance  of  o-  to  x-mode  coupling  near 
the  plasma  frequency,  which  may  result  in  intense 
radiowave  energy  absorption  in  the  ionospheric  plasma 
[41,46].  Theory  and  simulations  in  this  report 
indicate  that  mode  coupling  may  provide  a  viable 


explanation  for  many  experimental  observations  and 
should  be  an  efficient  heating  mechanism  for 
enhancement  of  future  experiments.  The  approach, 
conclusions  and  recommendations  regarding  mode  coupling 
and  its  relationship  to  high  power  HF  modification 
experiments  are  reviewed  here. 

In  Chapter  1  the  complex  nature  of  the  ionosphere 
due  to  the  Earth's  magnetic  field  and  the  importance  of 
high-power  HF  ionospheric  modification  and  its  possible 
application  to  over-the-horizon  communications  were 
briefly  discussed.  In  Chapter  2  the  relevant  theory  of 
electromagnetic  wave  propagation  in  the  ionosphere  w  ac¬ 
re  viewed  as  i  basis  for  understanding  tne  coupling 
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phenomenon.  This  chapter  considered  the  differential 
equations  governing  electromagnetic  wave  propagation, 
the' full  wave  solution  technique  and  the  method  of 
extracting  wave  components  from  the  total  field  to 
determine  the  reflection  coefficients.  In  Chapter  3 
parameters  of  the  three  ionospheric  models  chosen  for 
study,  demonstration  of  the  theory  with  sample  full 
wave  solutions  and  roots  of  the  Booker  quartic  and 
reflection  coefficients  were  presented.  The  major 
results  are  briefly  summarized  below: 

For  the  vertical  incidence  case  standing  waves 
below  X  =  1  indicated  the  presence  of  strongly 
reflected  downgoing  waves.  The  reflection 
coefficients  for  the  sample  model  were 
R  (  E )  =  0.889  and  R(H)  =  0.889.  The  imaginary 
part  of  q  was  very  small  throughout  the 
propagation  path,  indicating  small  energy  loss 
this  is  also  supported  by  the  large  reflection 
coefficients.  In  this  case  vertical  incidence 
ionospheric  heating  is  seen  to  be  very 
inefficient  . 

For  propagation  at  the  critical  angle  of 
incidence  in  the  3-N  direction,  there  were  no 


standing  waves  below  Xr,  indicating  the  absence 
of  any  strongly  reflected  wave.  The  reflection 
coefficients  for  the  sample  model  (Model  1) 
were  R ( E )  =  6.4450  x  1  0  ~  3  and  R(  )  =  6.9295  x 
1 0 —  3 .  The  roots  of  the  quartic  indicated  that 
the  upgoing  o-wave  couples  to  the  upgoing 
x-wave  which  then  reflects  just  above  the  X  =  1 
level.  This  x-wave  travels  downwards, 
approaching  the  plasma  resonance  region  Xr. 

The  imaginary  part  of  q  is  very  small 
throughout  the  propagation  path  except  in  a 
narrow  region  near  the  XR  level  for  the  x-wave. 
This  results  in  attenuation  of  wave  fields  and 
absorption  of  the  incident  energy  by  conversion 
to  heat.  Peaking  of  the  Ez  component  was 
observed  near  the  plasma  resonance  region. 

This  implies  that  the  absorption  of  energy, 
combined  with  the  large  Ez  field,  should  be 
capable  of  significantly  modifying  the 
ionosphere  near  the  Xp  level. 

For  propagation  at  the  critical  angle  of 
incidence  in  the  N-S  direction,  just  like  the 
5 - M  case,  there  were  no  standing  waves  cel: w 


the  XR  level,  indicating  absence  of  a  strongly 
reflected  wave.  The  reflection  coefficients 
for  the  sample  model  were  R(E)  =  2.4455  x  10-2 
and  R ( H )  =  2.3^33  x  10-2.  In  contrast  to  the 
S-N  results,  there  were  standing  waves  between 
the  Xr  and  1  +  Y  levels.  The  roots  of  the 

Booker  quartic  indicated  that  the  upgoing 
o-wave  couples  to  the  upgoing  x-wave  near  the  X 
»  1  level.  This  x-wave  then  reflects  from  the 
1  *  Y  level  and  approaches  the  plasma  resonance 

region  XR.  The  imaginary  part  of  q  was  very 
small  throughout  the  propagation  path,  except 
near  the  XR  level  which  results  in  attenuation 
of  the  fields  and  absorption  of  the  incident 
energy  through  conversion  to  heat.  This 
absorption  is  expected  to  modify  the  ionosphere 
significantly  near  the  XR  level. 

The  effective  coupling  region  around  the 
critical  angle  of  incidence  was  determined  for 
the  three  models  considered.  The  results 
revealed  that  the  coupling  window  is 
approximately  circular,  and  that  the 
determining  factor  affecting  its  angular  width 
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Experiments  need  to  be  performed  to  validate 
the  findings.  These  should  consist  of  ground- 
based  ionospheric  soundings  as  well  as  in  situ 
probing  of  the  electromagnetic  wave  fields  with 
rockets  and  satellites. 

Existing  experimental  data,  obtained  at  several 
locations  of  the  Earth,  should  be  examined  with 
the  emphasis  placed  on  the  o-  to  x-mode 
coupling  phenomenon. 

The  phenomenon  of  o-  to  x-wave  mode  coupling 
should  be  carefully  considered  in  the  context 
of  future  high-power  HF  radiowave  heating 
experiments  . 

The  theoretical  framework  of  the  coupling 
phenomenon  should  be  expanded  to  include 
nonlinear  effects,  such  as  plasma  wave 
generation  and  irregularity  scattering. 
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Reflection  Coefficients  for  the  Ionospheric 
Mode  Is  1  ,  2  and  3 . 
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Table  B.l  -  Reflection  Coefficients  -  Model  1,  in  the  Magnetic 
Meridian. 
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R .7  TORE ■ 
R .65*  l£* 
0.54  7?f  ■ 
7 .444  3E ■ 
6 .4  4  4  If ■ 
4.01 7?f- 
3.1 7436- 
1.37046- 
6.44406- 
I  .4  0  7?f ■ 

llYtFlf' 

7.3477E- 
8 .074  IE ■ 


HI 

l.R  077f  -01 
urr  i  lf-oi 
l.7fR3f-OI 
».*1  17F-0I 
I  .73  77  f  -0  I 
I  .6  1  481-01 
>.•0  716-01 
I. 77471-01 
•  6  6  ?0f  -01 
I  .f  9641-03 
I  .  774Rf -01 
>.470  11  -01 
.  71  07F  -01 
».*  7  4  4f-OI 
I  .  1  3  301  *00 


Ht 

R.Rf  77f  - 
0.77*7f  ■ 
0.4  76‘f- 
8 .4  ?4Cf  ■ 
7.07941  ■ 
6.3C4.76  ■ 
4.0440E ■ 
3 ,04 1 Rf ■ 
i .lisof • 
6 . 7  7  7  9f  ■ 
1.37  36F  ■ 
3.444  II- 
4.14106- 
7.I441F ■ 
7.670  IF  • 


Hi 

.00  77f- 

•  .050  11  ■ 
.0°63F ■ 

^ .90111- 
•.47941- 
» .06431 ■ 
1 . 3*  C6f - 

•  .34066 * 
I.47C4I- 
I  .9?94f ■ 
I  .5  107f  ■ 

t  .07301- 
f  .44131  • 


it  Cl 

)  70001  *00 
» .o.ccoi  *oo 
>  .open  *oo 
i  .non i *or 

1.00001 *oo 

i.-rcn  *co 

1.4C001 *oo 

i  .ocooi *01 
I  .04001  *01 
l  .07001  401 
1 .10001  *C1 
l  .1  6001  «Oj 
l  .70001  *C1 

i  .iccor  *c i 
i  .tcon  *o.i 


NAVE  0REOUENCT  •  5.00001*06 

HAG.  1  VECTOR  •  3.1590C-01 


GECHAC.  FIELD 


4.0470E-O5 


I-CIN  COS  Of  T  -  -3.70906-01 
T-01R  COS  OF  Y  •  0 .OOOOE *00 
/-BIN  COS  OF  Y  1  4./546E-0I 


CRITICAL  INC  ]  Cl  NCI 
START  INC  ME  ICHT 
ENDING  HE  ICHT 
INCREMENT 
SCALE  HE ICHT 
COLLISIONS 


-1 .07001 *01 
4. OOOOE  *01 
I  .COOOE  *03 
-4. COOOf -07 
I  .COOOE  *0* 
I  .OOOOE  *03 


REFLECT  I CM  COEFFICIENTS 


El 

0.88716-01  8 
8.9  T4Rf -0  1  8. 
4.074JE-01  8. 
4.103  IE— 0 1  0, 
0 • 79066 -0  I  0. 
7. 194  4C -0  1  6. 
*.63606-01  4 

l 

0  .7  386E  -a  3  7. 
I .64706 -01  I 
4 .04676-0 1  3  . 
8. 30446-01  7 
9.6  10  16-01  * 

1  .00446*00  S 


70716 
6404E 
4404  6 
0  34  f  E 
.44076 
0  4006 
1 0  34  6 
31496 
74696 
4  040  E 
460  IF 
1  3046 
04046 
11806 


6  T 

■01  0.88776-01 
■01  9. 01176-01 
•01  9. 1  4  706-01 
■01  9,7406f-0! 
-Ol  0  . 9744E-01 
■01  7  .3  I  4  46  -0  i 
-01  5.77316-01 
■01  3.67066-01 
•Ol  1. 4146E-0I 
-03  8.340  16-03 
■01  1.77446-01 
■01  4.19106-01 
•01  0.40476-01 
■01  0  .0  7716  -0  1 
■01  1.07896*00 


HI 

0.08776- 
8.70346- 
0.*?30f- 
0.’4O7F- 
7.71406- 
6.19676  ■ 
4  .00406- 
3.00106- 
1 .744*6- 
4.044  ]f- 
I .344  16  - 
1.3*006- 
6.6770F- 
7.473RF- 
7.35466- 


HT 

.08776-01 
.0*476-01 
.73106-01 
.4  It ’f -01 
.ICIE6-01 
•4EIC6-OI 
.07C6E-OI 
.  71706-01 
.4**96-01 
.67006 -03 
.77476-01 
.  3C  0  0F  -0 1 
.8ei*€-ol 
.075  7f *or 
.00616  *00 


HT 

F .80  77F- 
8.89046 ■ 
0  .0  *19  IF  ■ 
6.91  1 8  f  ■ 
0.44436- 
t. 86356 • 
4.3761E ■ 
?  .  30  76F  • 
1.41766- 
7.81176- 
1  .60196 • 
*.06046 ■ 
>.86096- 
0.90  !*f  ■ 
4.37436- 


inri 

.00006  *C0 
.00006  *00 
.OOOOf  *  no 

.aooof  »co 
.aoorf  *00 

.OCOCE  *C0 
.  4C00E ♦ CO 
.oron  *oi 

.04001 *01 
•  070CE  «Cl 
.  IOOOF  *0| 
.  I  *0CE *01 
. If COE *C  1 
.  4C00E *f  l 
.60006*01 


-  1  30- 


I 


WWW 


Table  B.2  -  Reflection  Coefficients  -  Model  1,  off 
the  Magnetic  Meridian. 


HAVE  fBEOUENCT  •  9.0000E*06 
M(.  r  V6CICB  •  3.I59CE-0I 
CECHAC.  MELD  •  9.69266-09 
I-OIR  COS  Of  T  •  -3.78906-01 
Y-OIR  COS  Of  Y  •  O.OOOOE  *00 
I-CIR  COS  Of  T  ■  9.29466-01 


CP  I T  1C  AL  INCICfNCf 
S  TAPI  I NC  HflCM 
FNC INC  HE  ICHT 
INCFEH6NI 
SCAIE  Mf  ICHT 
CClllSIPN? 


1  .C696f  *C1 
6 .COOOf *03 
l.COOOf *C? 

-A .COOOf - 0 7 
1  .COOOf ♦ C  A 
l.COOOf  *f  ’ 


6  .9^006  *0 
2.20081  -0 
4.23966-0 

l:UUhi 

8.946  U  -0 
8.69846-0 


RfflECTlCK  COfff  IClfMS 

fY  £1  M*  NY  HI  AHCl 

3  7.49906-03  6.27  796-03  6.4440f-03  7.8964F-03  6.9799f-f3  O.Or0Cf*<-n 

I  2. 49796-01  2.19896-01  2.6442 f -Cl  2.1013F-OI  2.3210f-Cl  9.00001-01 

1  4.8036f-01  4.19486-01  9.0840f-01  4.0444F-0I  4.4697r-CI  I  .00 001  *00 

i  i:mW\  S:n?*’i:$l  V.IWW I  UJWtf 

I  3.««6ir-OI  8.3767E-OI  9.<»44t6-0l  5. DBM  -Cl  f  .9  3  A  71  -C  1  9  .  nor  01  •  r  0 

1  9.!OUf-OI  8.48906-01  9.60971-01  8.3790F-01  f.9A**l-OI  B.nooo««ro 


WAVE  FBCOtlENCY  •  9.00006*06 


MAC.  T  Vf  C  TCP  • 
CECHAC.  f  IE  ID  • 
*-01R  COS  Of  T  • 
T-01P  COS  OF  T  . 
Z-flIB  COS  OF  r  • 


•  3.1990E-OI 


■  9.64266-09 

•  -3.7890E-01 


O.OOOOE  *00 
9.29466-01 


CRITICAL  INCTCfNCE 
STAPTING  ME ICHT 
f NO  INC  HEIGHT 

ihcpement 
SCALE  HEIGHT 
CCLLISIONS 


-I.C700E *01 
«.0000f  *C3 
l.COOOf  »03 
-9 .OOOOf -0? 
I.OOOOf *0A 
I  .COOOf  *03 


El 

8.23866-0 
2.64246-0 
A  .80926-0 
7.91806-0 
4 .1 90  46 -0 
9.39196-0 


PEflECTICN  COEFfICIENIS 


9696-03  8.69766-03  8.3981f-03  6.89916-03  7. 8117F-03  O.OCCOFMO 
.  . .  2.199AF-0I  7.767*6-0!  2.9C566-0!  S.OCOOE-Cl 


29696-03  8.69766-03 
.  -.32826-01  7.69366-01 
I  A.7A076-01  A. 90766-01 
I  7.00A7E-01  8.07316-01 
1  8.1  3616-01  9.33  I  76-01 
1  8.48466*01  9.98896-01 


A • 006  9f -0  I  9  .036‘E-OI  A.96lAf-0| 
6.623  lf-01  8.790A6-0I  7.9IP4F-0I 
7 ,70  3  If-01  9.97686-01  f.7C»Af-Cl 


6  .069  96 -0  I  9.82076-01  8.989BE-0I 


9  6  I  Af -0  I  l.OCOOfACO 
9  I 8Af -0  I  2 .00001 ♦ CO 
7CA4F-C1  3.0C0Cf»''0 
98986-01  9  .OOOOf  *00 


9.2BS1E-0I  8.69 366-0 1  9.A702E-OI  8.2976f-0l  9.69226-01  f.9P93f-Cl  8.0C0C6M0 


-  1  31- 


«  r.' 


iV 
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Table  B.3  -  Reflection  Coefficients  -  Model  2  -  in  the 
Magnetic  Meridian. 


HAVE  FREQUENCY 
MAC.  Y  VECTOR 


5.00006*06 

3.15906-01 


CECRAC.FIELO  •  5. 69266-08 

X-CIR  COS  OF  Y  *  -3.7890E-01 

r-oia  cos  of  r  •  o.ooooe*oo 

2-01*  COS  OF  T  •  9.75966-01 


CRITICAL  INCICENCE 
STARTING  H61G6-T 
6  NO  I NC  ME IGHT 
INCREMENT 
SCALE  HEIGHT 
CCLL IS  IONS 


1  .06-666*01 
]  .COOOE  *0  3 

o.ooooe  *r o 
-i  .ooooe-oi 
I  .COOOE  *03 
I  .COOOE *03 


REFLECTION  COEFFICIENTS 


EX  6  T 

9  .22596-01  9.27*96- 
6.99986-01  9.1097E- 
3.16916-01  8.9*956- 
8.63116-0!  7.0OO5E- 
A .29  22E-0  I  A  .56196- 
3.52366-01  3.80976- 
2.7659E-0I  3.009*6- 

i:?2ill-ci  f:J;A?f: 

A. 35696-07  9.I917E- 
•.19876-02  9.11056- 
7.61ACE-0?  7.8875F- 
7 .299  IE -0?  7.59766- 
3.07AIC-02  3.9010F- 
*.901*6-0?  6.5*796- 

JrfcWf: 

3 .700  IE  *0  I  9  ,»99|6- 
8. IT7T6-OI  7.MM6- 
7.78C3E-OI  9.11776- 


ET 

■01  9  .2 2S Af  • 
•01  8.99106' 
•01  8.16996 
•01  6.69016 
■Ol  9  .29  79F 
•Of-  3 .92876 
■01  7.76966 
•01  1.97726 

-01  I  .  1 6  106 

•07  M. 37066 
-07  *.20776 
-O?  7.61866 
-07  ?. 30396 
-02  1.0  1966 
•07  5.917*6 


01  1  .90006 
Ol  7. 20996 
01  1.70A16 
01  6.18976 
01  7.80966 


HI 

■<71  9.229  96-01 
01  9 .169  IF  -  OL 
■01  8. «6S1E-0I 
•0  7.l’3«F-OI 

•0  9  •#  4  ??F  -0  i 

■0  1.91936-01 

■0  3.09386-01 

:§  ?:«?5WI 

■02  9.80776-07 
■C2  9.92716-02 
■07  ?.99]Af-07 
■07  7. 63*76 -0? 
■07  1. *7876-07 
•07  6.79176-0? 
■01  I.67066-0I 
■ol  J.*7I?F-Ct 
•01  9.1919f-0l 
■01  7. *8716  -01 
■01  9.99796-01 


HT 

9. 72596' 
8.90116- 
8.07976' 
6 .52976 • 
9.19896- 
1.A9C76 ■ 
2.69366- 
1.92I96- 
1.17876- 
8.1 181E' 
5.09  796' 
7.5’6’F' 
7.23016' 
2.9FC96 
5.71196' 
1.1*706' 
7. 178  ’6- 
1.8M  IF- 
5.909*6 ■ 
7  •  1887F  ' 


HT 

■01  9.77596 
01  8.07916 
01  8.30C86 
0!  F. 80096 
01  9.38*96 
01  3 .65026 
01  2.87076 
Cl  2. 05166 
ol  I.21C36 
07  8.71816 
0?  5.97616 
0?  2.73176 
07  J.AC’96 
C7  ?.?l«?f 
0 7  8.1 7796 
01  1.96‘Nl 
01  2.1I76F 
01  3.906IF 
01  t.  5  7716 
Cl  8.35*06 


■Q  I  0. 
0  1  7. 


oi  a. 

•01  9. 
■Ol  9. 

■ol  1  . 
•0?  1. 
•02  l. 
■C  7  1  » 
■0?  I  . 

c?  I. 

C?  I. 

ol  I: 

■ol  1. 

■Ol  I  . 

•Ol  l. 


ANCt 

00006  *00 
OCOCF  *00 
00006  *fO 
crrrf  *rr> 
OOOCE  *00 
5C0C6  *C0 
00006*00 
*cecF*co 

0C006  *0 1 
07006  *C» 
09006  *01 
0FCC6  *CI 
0700F  *01 
08C0F  *  r  | 
iocor*oi 
j  *006  *01 
?ficct  *r i 
lrcCE  *c| 
*  CO  OF  *01 
7  CODE  •  C  l 


HAVE  FREQUENCY  «  5.0000E*06 

N8C.  »  VECTOR  -  3. 15906-0 1 

CECMAC.  FIELD  •  5.6926E-0* 

1- 018  COS  OF  T  •  -3. 78906-01 

Y-0I8  COS  OF  Y  •  0.00006*00 

2- 018  COS  OF  Y  •  9.7596E-01 


CRITICAL  INC  I CENCE  •  -1  .07006*01 


START  INC  HEIGHT 
ENDING  HEIGHT 
INCREHf NT 
SCALE  H6ICHT 
COIL ISIPNS 


1.00006*03 
C. 00006  *00 
•1 .COOOE -Cl 
1 .COOOE *03 
I .cooor  *03 


8EFLECTICN  COFFF ICIENTS 


61 

9.22996- 
9.10006- 
8.A9J3I- 
8 .9  7596 - 
9.59776- 
3.80096  - 

hKUK 

*  .96706- 
3.01 I7f- 
7.99996- 
i. 12926- 
F.  17056- 
1  .51396- 
?.AI8Cf- 
9.19336- 
7.0989F  - 
9.1  3996  - 


6  Y 

01  9  .77596- 
01  8.99086- 
01  8.19666- 
01  6.60966- 
Ol  9.77766- 
01  3. 51556- 
01  7 . 7f 926 • 
01  1.98776- 
«l  1.17836- 
0?  8 . 5*6?F- 
07  5 . 9 C 88  6 ' 
0?  2.77696- 
0?  2.7  1*56 ■ 
0?  7 .8/676 ■ 
0?  *>.56  706' 
01  1.15676' 
0)2. 1*566 
01  3.65176' 
0  16. 098*6' 
01  7.61776' 


01  9  . 
01  9. 


01  9  . 
01  3. 


01  9  . 

01  7  . 


6  T 

22596 ■ 
0968F- 
9 159F' 
96836 ■ 
59166- 
79996- 
999  76- 
I  6  I  76  • 
79126  • 
39566- 
95  186' 
00696' 
99556- 
1?  366' 
18916' 
510*6' 
91  176 
1  1996 
0  7*76  ' 
10696 ' 


M 

01  9.775«F-01 
01  8.8860F-C1 
01  8  . 09676-01 
01  6.98806-01 
Oi  9  .  |  70  3F  -C  I 
01  3.9707T-C1 
01  7.F8*3r-fl 
Ol  1.92  1 9f  -«  I 
01  I.IAOTF-OI 
0?  8.77186-02 
0?  * .779CF-C? 
0?  7  .f  3  1  76  -0? 
0?  7. 137 76 -G? 
07  2  .72691  -C2 
07  5. 36096-02 
01  1 . 30  9  9f  -f  1 

0}  ?.P6*«f-CI 
01  1.98916-f 1 
0  1  5 .77901  -0  1 
01  7.I98F6-CI 


HY 

9.77*96-01 

9.19996-01 

8  .*  3996  -01 
7.09776-01 

9  .6*0*6-01 
3  .89706-01 
3.081  7F-01 
7.72196-01 
1  .130*6-01 
9.68786-0? 
6.1 9C*f -0? 
3.10906  *0? 
7.57*76-0? 
3.727C6-0? 
6.36796-07 
1  .56976-01 
7. 80306-01 
9 . 1C  796 -01 
7  .9  1  “F  -C  1 
« .6  I  566 -0  I 


HI 

5.77*96  ■ 
9.07996- 
6  . 1 C 79  6  - 
t  ,8f  *06- 
9.35926- 
3.66866- 
7.85906- 
7.06191  - 

1 .2907r- 
5.07106- 
* .7f59r. 
7. 88  196- 
2.39  ’36- 
}.99iar- 
5.69976- 
1 .99796- 
7.3CC36 • 
3 .97* OF- 
I  .K  ‘66- 
6.5  0  ?06  ■ 


if  ri 

l.OCOOf  ♦«  0 
>.00006  «co 

>  .00006  *C(> 

>  .00001  *co 
I.oooo6*ro 

1.8COOI  *00 
6.00006  *6  0 
*.‘0006  *<0 
1  .C0006  *n 
I  .07006  *01 

1 . 09  co6 *r  i 
[  .of  on  *r  i 

I  .07006  «C  I 
l  .06001  *r i 
1  . I  0006  «f I 
[  . | ‘006  *01 
I  . 7  C  006  «  C  I 

1  .  ’0006  *0  I 

I  . 5000*  ♦  r  i 
I  .  70006  *0  1 


'WZ* 


W3 


V*V-  *>  /i  i  v  .■  V  V '/.V 


^v,mww.v  ^  u  VW  MU  J  JW  i^gyu  v^r^rz  a'^.v.-v? 


Table  B.4  -  Reflection  Coefficients  -  Model  2,  off  the 
Magnetic  Meridian. 


NAG.  T  VECTOR 
CECHAC.  FIELD 
I-B1R  COS  OF  T 
Y-BIR  COS  OF  T 
I-BIR  COS  OF  T 


3.  I940E -0  I 
9.6426E-09 
-3.789CE  *0  1 
C.OOOOf *00 
4.2946E-01 


STARTING  HEIGH 
6  NO INC  HEIGHT 
INCREMENT 
SOLE  HE  I  CHI 
Cfll IS1CNS 


I. 00006*03 
0  •  OOOOf *00 
-1  .COOOf -01 
1  .  COOOF  *Ci  3 
J  .COOOf  ♦O’ 


REFLECTION  COEFFICIENTS 


4  .229  IE 
6  .OOt  OE 

2.009CE 

3.47111 

ml 


FT 

07  7.7167F 
07  3.87466 
07  9.43176 

8i  urn 

at  1.8146E 
01  3.J177F 

si  }-l»H 


FT 

0?  7.44636-07 
07  4.2717E-0? 
€7  9.14496-07 


81  UlttlW] 

01  7.00176-01 
I  3 . 661 4F -01 
1  '.1171F-0I 

I  S:i;iU:gI 


FI 

7.I3F4F-0? 
3.6904E-0? 
1  .74  I  Of  -0? 
4 .47707  *07 
I •  74  41f  -0  I 
I .71177 -Cl 

3  .71731  -01 

4  .■>1411  -01 


H7 

7.1764F 
4.3600E 
A ■ 14 |4F 

\:mi 

7.0664E 
1.77#  If 
i.iinr 

8:8ISC4 


HI 

0?  ;.34»Of 
07  4 • 04 “ If 
07  *.7440f 
-0!  i.c#?ir 

-01  1.4  I  f  9f 

-01  1 .4  I P4  f 
3.4  I  Ilf 

4.4f>4f>l 
#.77041 
E  .01 13f 


-01 

-cl 

01 


Al.fl 

-o?  o.ooo of *ro 

-07  I. 0000# -Cl 
-07  2 .00001 -Cl 

--?}  *:S8S8f:n 

-Cl  1  .00001  *00 
l.oroni *ro 
l.oocoi  *oo 
4 . ocoot • ro 

6.00001 *C0 


:$! 

:?! 


HUE  FREOUENCY  •  1.00006*06 
MAC.  Y  YE  C TOO  •  3.1410E-01 
CECHAC.  FIELD  >  1.64  76E-0S 
I-eiR  COS  OF  T  •  -3.7040E-01 
Y-Mt  COS  OF  Y  •  0.00006*00 
7-BIR  COS  OF  T  •  9.2444E-0 1 


CF II  1C  41  INCICENCE 
STAR  7 1 NC  HEIGHT 
6 NO  INC  HEIGHT 
1HCREHENT 
SCALE  HFICNT 
CriLISIONS 


1  . C  6466 ♦ C  1 
1  .  COOOE *0  3 
0  .COOOE  *00 
-I. COOOF -01 
i  .rooof  *oi 
I  .00006*03 


REFLECTION  COEFFICIENTS 


61 

7.711  IE 
7.1666F 
I  .01616 

1.38136 

7.44646 

t'.im 

7.4IIEE 


f  T 

07  7.91616-07 
03  7.17436-03 
0?  I .17116-07 
07  6.46016-07 

02  l.aoiof-o] 
01  I .13906-01 
01  3. 29076-01 

81 2:fWf:S! 

at  a .7 199E-01 


67 

7.79716-07 
7.1 743E-03 
1 .01116-07 

9 Io8o^l -82 

1 .34  166-01 
2.91296-01 
4 .394C6-0J 
9  .66616-0  I 
7.60836-01 


HI 

2.6289F 
8.21296 
I  .16416 
6.69166 
I  .03446-' 
1  .99066  - 
3.3661F- 
4  .48116  - 
6  .3116#- 
8  .46606- 


HT 

02  2.22496-02 
03  6.991E6-03 
02  9.89496-03 
02  4.66906-07 
1  1.80766-07 
l  !. 34776-01 
01  7.86126-01 

81  ?:lo6UI:81 

01  7.41076-01 


N  1 

7. 39106-07 
7.49236-01 
1.06776-02 
t .10136 -C  7 
9.41*66-07 
I. 41206-C1 
3.0  1936-01 

4.17736 -Cl 
«.844«F-f  I 
7.88616-01 


ANCL 

.OOOOf  *ro 

.ocooe-ci 
,00006 -0 1 
,  OOOOf -Cl 
.OOOCt -r 1 

.ooooe  *ro 
.  OCOOE ♦ CO 
•gccci *oo 
.OOCOf  *ro 
.00006*00 
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Table  B.5  -  Reflection  Coefficients  -  Model  3,  in  the 
Magnetic  Meridian. 


WAVE  FR60U6NCT 

mac.  r  vc cron 


8.00006*06 
3.1 590? -0 1 


G60HAC  •  F  ICtO  •  9.64266-08 
l-om  COS  OF  V  •  -3.78906-01 
V- 01*  COS  OF  T  «  0.0000E *00 

i-em  cos  of  v  «  9.29466-01 


critical  iNcirFuce 

STARTING  HF IGF  T 
FNOING  HE 1GFT 
INCREMENT 
SCALE  HEIGHT 
COLLISIONS 


I .C696E  *C  1 
I  .COOOE *03 
C .OOOOE *00 
-2.C0006-0? 

1 . COOOE *03 

2 . COOOE *07 


REFLECTION  COEFFICIENTS 


El 

9.H69F- 
•,0379*- 
8.24406- 
6.7000E- 
4  .29886- 
3.97306- 
2.80946- 
2.01C6E - 
1.  19  6  76- 
8.99186- 
*  .107  *6* 
2.812C6- 
2.09166* 
2  .82  766  - 
* .76  I  IE” 

i  :8«|: 

3.79716- 

6  .7917F  - 

T.asoif- 


01  8. 

01  7. 


01  3. 
01  7. 
01  I  . 
•07  9. 
02  S. 
02  2. 


FT 

3I69F- 
l»Mf* 
947TE- 
0709F- 
C7tnf - 
8676F- 
08776- 
19866- 
10766- 
40676- 
8446E- 
7749F- 
7667F* 
I  3016- 
*91 ?F- 
S444E - 
49  74  F  - 
788<7f  - 


ET 

.31696* 
.04  |0E 1 
.2497F* 
.  70696* 
.  30  166* 
.81876* 
.81  366  * 
.01 14E" 
.18876- 
•86606 1 
.31 I4f • 
•8TC4F* 
.099  IE  ■ 
•  8  7  2  IF  * 
•77186* 
•40826- 
.2261E- 
.7604F 


01  9. 
Ol  9. 
01  8. 
0!  7. 
014. 
01  3. 
01  3. 
•01  7. 
■01  I  . 


02  3. 
02  6  . 

ol  ): 

01  4. 
01  7. 
Ol  I  .. 


HI 

3164E-01 
78896-01 
648  7F-01 
709  IF -0  I 
74  |  36-01 
O707E-0I 
1430F-0I 
7780E-01 
346*E-ol 
729«E-02 
0*0  I  F-0? 
8746 f -0? 
38I7F-P? 
7481F-C? 
631 7F-C2 
6266F-C1 
9*67f-0i 
48941-01 
6"06»  -f  | 
001  If  *00 


HT 

9.3169F- 
•  .98866  - 
8.18  4CF ■ 
6 .86  4°F  * 
4.7CC8f • 
3.4689E- 
7.7397E- 
1.9*74F- 

I.1976E- 
8.1076E- 
8.1479F- 
2.4  374F- 
1.99CCE- 
2.741<F- 
9.883CF- 

I. 387CF- 

J. 1464E- 
1.6  17  If  ■ 
«  .OF67F  - 
/ .4 -80E ■ 


I.OOOQF *CQ 

l.oooof  «r.A 

>  .orooE  *ro 
i.OOOCf  *co 
I.OCOOf  *ro 
I.8C0CE *C0 
> .OOOCE  *C0 

i.srooE  *ro 
[  .ncccF  »c i 

I  . 070PE *0 | 
l  .0400F *f I 
[  .06001  *ri 
l  .07C01 *01 
l  .OPf.Ol  *0  1 
l . lOCOf  *01 
I  .1*001*01 
[  ,2r6oi *c I 
i  .icon  *m 
l  .scori  •'•1 

I  .  7 C001  *C| 


HAVE  FREQUENCY 
MAC.  T  VI  CTO* 


9. OOOOE *06 
3. 18906-01 


CFITJCAI  INCICFNCI  •  -1-07006*01 

STARIINC  HEIGHT  *  1.00006*03 


CECHAC.  rlELO  ■  9 .64266-08 
8-OIR  COS  OF  T  •  -3.7890E-01 
V-C1R  COS  OF  T  •  0.00006*00 
2-CIR  COS  OF  T  •  9. 29466-01 


STAR  1 1  NC  HEIGHT 
f HP  INC  HEIGHT 
TNCFCHENT 
SCALE  HC1CHT 
CrtllSIPNS 


0.00006*00 
-2. COOOE -07 
I  .COOOE *03 
2 .COOOE  *02 


REFLECT ICK  corn ICJf NTS 


El 

* ,06  3CE " 
9.41CCE- 
8. 41746* 
7.01846* 
4.61 346* 
3, 86066* 
3.0888E- 
2.70496* 
1  .3  7  0  36* 
8.62896* 
6 .179  IE* 
3  o  108  4F  • 
7.4803?* 
-  3. 12  4  66 ■ 
6.17706* 
1.82416* 
7  .4  38  36  * 
■  4.179  l£* 
7 . 1  87  IE  ■ 
9  .  2  1  0  76* 


1  T 

01  9. 00306-01 
01  8.88746-01 
01  8 •  12  36 F -01 
01  I .69076-01 
Of  4 .7886F-01 
01  3. 87126-01 
01  2.81746-01 
01  2.01926-01 
01  1 .20286-01 
02  8.79386-02 
07  9.98736-07 
02  7.81446-02 
0?  7.29026-0? 
02  2.82766-02 
07  8 .8679F-07 
01  1. 36636-01 
01  2. 1737E-0I 
Of  3*6832  £ -01 
01  6. 1810F-01 
01  7. 69816-01 


f  T 

9.0630F- 

9.00696- 

6.41196* 

Miffi: 

3.899  1E- 
3.09  106- 
7.20146- 
1.31816- 
9.6122E- 
t.l 191E- 
3.10*96- 
7. 48*96- 
3.11906- 
E  .  16886- 
I.8717F- 
7.433SC- 
-4.1  7C7E- 
7.13*86- 
9.18206- 


F« 

.01  3CF  ■ 
.  79  8  7f  ■ 
•0  2  3  4f  ■ 

:)V4l 

.  4  7  44  F  * 
. 73 l?f • 
.9*771* 
.11 3PF ■ 
.46711 • 
•  3t  0* F  * 
.7167F* 
.171 76 • 
.72241* 
.  It  6CF  ■ 
.  3  1 4  ?F  ■ 
.006  3F  • 
.81996' 
.87  34f  • 
.7C77F 


Ot  4. 
0!  3. 
01  3. 


07  9. 
0?  I. 
07  3. 
07  7. 
02  3. 
07  6  . 


FT 

01  30F  -01 
OS  94F  *0  I 
81C3C-01 
I  38  PE  -0 1 
71 76f -Oi 
9833C  *01 
I  339F -0  1 
71 4*f -01 
388  ?F -Cl 
91176-02 
3C94f  -0? 
7C4  1F-C? 
868  3f  -07 
222  36-07 
374CF-07 
8  78  7f  -0 1 
8  74  IF -0 1 
34076-01 
4830F -01 
69896 -0  I 


Fl 

‘-.Of’Of 
E.9386E 
I • 2  7F°F 
I .84-9F 
4 .46776 
3 . 7  ?<  8 F 
7.94  ?8F 
7.1 1 96f 
1.76666 
S. 2  7  SOf 
‘.81 I6F 
7 .9  749f 
7.38C0F 
2.98786 
8.4009F 
I  .4  8  30F 
7.3196F 
3.9894E 
1.7 ??6f 
1  .8  771F 


01  0. 
CI-7. 
01-4  . 
01-6. 
0  I-fl  .. 
C  1-B. 
0  1-9., 

of-9. 
01-1  . 
02-1  . 
07-1  . 
C  7-1  . 
07-1  . 
02-1. 
07-1  . 
C  I  -  l  . 
01-1  . 
■01-1. 
01-1  . 
C  1-1  . 


a  nr  i 

proof  *or, 

OCOOt  *00 
OPO**f  *00 

ooooi *on 

OCOCl *00 
8C0CI *C0 
OCOOt  *rn 
80COI «C0 
OOOCf  *01 
02CCI *C1 
04001 *01 
0  70  Of  *01 
0  7001  *01 
Of  OCf  ♦  C  I 
ICCOf  *C1 
1*007  «C1 

2  ooor *ci 

3C0CE *C1 
8000f  *01 
70001  *C  1 
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Table  B.6  -  Reflection  Coefficients 
Magnetic  Meridian. 


Model  3,  off  the 


HAVE  FREQUENCY 

mac.  t  vector 

C6CHAC.  F 1610 
1*011  COS  OF  Y 

y-oir  cos  of  r 
j-cir  cos  OF  T 


5.00006*06 
J. 15106-01 
5. 64266-0S 
-3.76106-01 
0.00006*00 
1.25466 -0 1 


CRITICAL  INC  I C6NCE 
STARTING  M6 TCFT 
ENDING  HE IGHT 
INCREMENT 
SCALE  HEIGHT 
CULTS  IONS 


I .06466*01 
I  .00006*03 
O.COOOE  *C0 
-2.00006-02 
1.00006*03 
2 .00006*02 


REFLECT ICN  COEFFICIENTS 


El 

2.04716 
4.3831* 
1  .24606 

hm 

\'Axn 

4.44416 


FT 

•0?  2.26406-02 
•03  4.64716-03 
•0|  1.374(6-02 


-02  6.74446-02 


•02  1.034(6- 
•8l  1.57616- 


:S! 


°1 


1.57616-01 

3.10416-01 

4.(0626-01 

7.40106-01 

4.42446-01 


E  T 

.04  Ilf-02 
.34146-03 
,24026-07 


-07  6.14766-02 
-01  4.42316-02 
-01  1. 42416-01 


.00216-01 
.4  4  3C6-0 1 
.81  306-01 
.85546-01 


HI 

7 .34  6  IF  -0  2 
5.02341-03 
I  .43006-02 
7.0*726 -02 
1 .07746-01 
1  .63346-01 
3.47316-01 
4  .056  7F  -01 
7.64146-01 
4.66406-01 


HI 

1 .4(616 
4.25726 
I  .21046 
4.46766 
4.12576 
1.364CE 
2.404  16 
4.11836 
6.62446 
8.6634E 


HI 


■07  7 . 1  4C66  -0  2  0 
03  4.5e286-03  1 


*6-02  1.3C46E-C2  2 


-07  e. 47476-02 
-02  4.(3316-07 


5. 


•01  I. 44116-01  1 
•01  3.1  3C4f -0 1  2 
■01  4.63426-01  1 
■01  7.07676-01  5 
•01  “.12146-01  a 


ANCl 

00006  *C0 
QCOCE -0| 
OCOOE -01 
OCOOE -0  1 
OCOOE -C  I 
000CE *00 
00006*00 
occoe  *ro 

00006  *C0 
00006 ♦ CO 


HAVE  FREQUENCY  ■  5.00006*06 

HAG.  I  VECTOR  •  3.15406-01 

GECHAC.  FIELD  •  5.64266-05 

1- CIR  COS  OF  T  •  -3.78406-01 

r-CIR  COS  OF  T  *  0.00006*00 

2- 018  COS  OF  V  •  4.25466-01 


CRIT1CAI  INCItlNCI 
STARTING  HI  IGHT 
ENDING  HE  I  CHI 
INCREMENT 
SCALE  HE  ICHT 
CCLL I S IONS 


•1.07006*01 
1  . C 0006  *0  3 
0.00006*00 
*1 .COOOf -01 
1  .  COOOi *01 
2.00006*0? 


REFLECTION  COEFFICIENT! 


61 

7.48746 
4.28316 
6.076  36 

!. 14366 
.44846 
.025  16 
3.70846 
5.23146 
7  .6 1 4  46 
4.44706 


FT 

>07  2.2*446-02 
■07  3.(7346-02 
•07  5.44446-02 
■01  1.03446-01 
■01  I.3456E-01 
-01  1.837(6-01 
•01  3. 16446-01 
•01  4.76126-01 
•01  7.01576-01 
■01  6.14746-01 


6  T 

7.4(316-0? 
4.27566-02 
f .06EC6-0? 
I .14 I 76-01 
415(6-01 


1  .415(6-01 
7.02116-01 
3.70386-01 
5.22616-01 
7.623C6-0! 
4.57276-01 


HI 

7.170*6-07 
3.73  7  4  6  -  02 
5.30266 -07 
4.68166-07 
1  .30(16-01 
I  .7t8 76-01 
3.74836-01 
4  .60246  -Cl 
6.76576-01 
8.72146-01 


HY 

7.5E  446' 
4.4  1*76  • 
6.264EF 
I .1 7606 
1  .54476 
2.08766 
3.82266 
5.3(446 • 
7.840  IF  1 
4. 73436' 


Hi 

07  7.37436 
C?  4.04666 
0?  •  .81746 
01  1.04406 

‘  1.4  3341 

1  .43766 
3.55  176 
5  .01(46 
7.338BE 
4.23516 


AN6  I 

-c ?  o.oroof*ro 
-02  1.00006-01 
-02  7 . OOOOf -c 1 
-0  1  5  .00001  -01 
-Cl  7.0C00I-O1 

-ci  l.ooooi  *rn 
-ci  ?.ococf*'-o 

-01  3.00001*00 
-Cl  5.00006 *ro 
-01  8.00006*00 


E 

E  0 


